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Heat  Transfer  From  Water  Flowing  Through 
a  Chilled-Bed  Open  Channel 


PAUL  W.  RICHMOND  AND  VIRGIL  J.  LUNARDIN1 


INTRODUCTION 

The  determination  of  heat  transfer  coefficients  for  water  flowing  over  ice  has  been  of  interest  for 
a  number  of  years.  Initially,  research  in  this  area  developed  from  curiosity  as  to  the  cause  of  ripples 
on  the  underside  of  ice  covers  on  frozen  rivers  and  from  the  desire  to  predict  ice  cover  growth  and 
decay.  Recently,  work  in  this  area  has  been  conducted  to  obtain  heat  transfer  correlations  applicable 
to  the  use  of  large  masses  of  ice  for  dissipating  waste  heat  in  self-contained,  underground  installations 
(Lunardini  et  al.  1986). 

Observations  have  shown  that  the  heat  transfer  is  much  greater  for  water  flowing  over  ice  than  for 
water  flowing  over  flat  plates  without  melting.  Early  work  in  this  area  consisted  of  observations  of  the 
occurrence  and  behavior  of  wave  shapes  on  the  underside  of  river  ice  covers  (Carey  1 966, 1 967,  Larsen 
1969,  Ashton  and  Kennedy  1970, 1972).  Following  this  work,  a  number  of  investigators  (Gilpin  et  al. 
1980,  Ashton  1972,  and  Hsu  et  al.  1979)  examined  how  these  ripples  are  formed  and  their  effect  on 
the  heat  transfer  coefficient.  In  these  studies,  heat  transfer  rates  from  initially  plane  ice  surfaces  were 
given  limited  analysis,  the  greater  interest  being  in  the  formation  of  rippled  ice.  Hirata  et  al.  (1979a 
and  1979b)  examined  steady-state  melting  of  ice  on  plates  in  the  entrance  region,  and  confirmed  that 
higher  heat  transfer  rates  occur  than  with  plates  without  ice.  In  their  turbulent  regime  work,  where  the 
concern  was  with  the  mechanism  of  transition,  the  heat  transfer  rate  was  also  higher  than  predicted  by 
flat  plate  theory.  Haynes  and  Ashton  (1979)  presented  new  data  for  water  flow  in  large  aspect  channels 
and  compared  their  data  with  one  analytical  and  one  empirical  model,  without  good  agreement.  Zisson 
(1984)  conducted  experiments  in  a  large  open-channel  flume,  studying  the  melting  of  horizontal 
bottom  ice.  These  experiments  were  of  short  duration  (7-8  hours)  and  average  heat  transfer  rates  were 
determined  throughout  the  flow  regime,  including  some  natural  convection  experiments.  Further 
analysis  (Lunardini  et  al.  1986)  of  the  data  showed  that  for  low  flow  rates  and  with  the  water 
temperatures  above  3.4°C  the  heat  transfer  rate  does  not  drop  below  the  value  for  free  convection 
(approximately  489  W/m2).  In  general,  they  attributed  the  heat  transfer  rates  greater  than  the 
nonmelting  case  to  increased  free  stream  turbulence,  thermal  instability  due  to  the  density  inversion 
of  water  near  4°C  and  turbulent  eddies  associated  with  the  generation  of  a  wavy  ice  surface  during 
melting,  especially  at  high  Reynolds  numbers. 

The  present  study  was  initiated  in  an  attempt  to  further  define  the  cause  of  increased  heat  transfer 
in  waterflowing  over  ice.  The  approach  taken  in  this  study  was  toeliminate  some  of  the  characteristics 
of  melting  ice  such  as  turbulent  eddy  generation  due  to  a  wavy  ice  surface  and  flux  of  meltwater,  while 
maintaining  a  similar  temperature  range.  In  order  to  accomplish  this,  a  small  open-channel  flume  was 
constructed  for  which  the  heat  flow  from  the  bed  without  the  presence  of  ice  could  be  measured.  From 
these  data  it  was  anticipated  that  a  more  complete  understanding  of  the  heat  transfer  could  be  obtained. 


BACKGROUND 


Forced  convection  heat  transfercorrelations  are  generally  presented  as  a  relation  between  Reynolds 
number,  Nusselt  number  and  Prandtl  number. 

The  Reynolds  number,  a  ratio  of  viscous  and  inertia  forces,  is  used  tocategorize  flow.  The  definition 
of  the  Reynolds  number  is 

Re  =  ¥L.  (1) 

vr 

where  Re  -  Reynolds  number 
V  =  mean  fluid  velocity 

vf  =  kinematic  viscosity  of  the  fluid  at  the  film  or  bulk  temperature 
L  =  characteristic  length. 


Two  characteristic  lengths  can  be  used  in  correlating  heat  transfer  data  for  the  horizontal  flat 
surfaces  of  a  flume.  These  are  the  length  from  the  flow  inlet  to  the  point  of  interest,  and  the  hydraulic 
diameter,  defined  as  four  times  the  flow  area  divided  by  the  wetted  perimeter.  This  definition  of 
hydraulic  diameter  is  standard  for  heat  transfer  correlations  and  is  four  times  the  value  often  used  in 
hydraulics. 

The  Nusselt  number  is  defined  as 


Nu  =  hL  (2) 

kf 

where  h  is  the  convective  heat  transfer  coefficient  and  it  is  the  thermal  conductivity  of  the  fluid. 

The  Prandtl  number  is  the  ratio  of  the  diffusion  of  momentum  to  the  diffusion  of  energy,  in  a  fluid, 
and  characterizes  the  relative  growths  of  the  momentum  and  thermal  boundary  layers.  The  definition 


is 


(3) 


where  vfand  afare  the  kinematic  viscosity  and  thermal  diffusivity,  respectively. 

Although  these  terms  were  defined  by  material  properties  evaluated  at  the  film  temperature,  they 
can  also  be  defined  using  properties  based  on  the  bulk  flow  temperature.  Bulk  flow  properties  are 
generally  used  in  specifying  internal  flow  while  film  temperatures  are  typically  used  when  describing 
external  flows.  This  convention  is  used  in  later  analyses. 

The  data  obtained  in  these  experiments  were  for  the  flow  of  water  in  an  open  channel  with  turbulent 
flow.  The  Reynolds  numbers  (based  on  hydraulic  diameter)  extended  from  fully  turbulent  values 
through  Reynolds  numbers  generally  described  as  the  transition  region  values.  Bulk  water  tempera¬ 
tures  ranged  from  5°  to  33°C  and  channel  bed  surface  temperatures  ranged  from  near  zero  to  about  9°C. 
In  order  to  determine  the  effect  of  the  density  inversion  on  the  heat  transfer  coefficient,  data  were 
collected  for  bed  surface  temperatures  below  4°C  and  were  compared  with  data  collected  at  bed 
temperatures  above  4°C.  It  was  expected  that  the  data  obtained  with  the  bed  temperature  below  4°C 
would  compare  well  with  those  for  flow  over  melting  ice.  Similarly,  the  data  obtained  with  higher  bed 
temperatures  should  compare  well  with  standard  heat  transfercorrelations.  Heat  transfercorrelations 
for  melting  ice  and  for  water  at  low  temperatures  are  limited.  The  majority  of  the  correlations  for  water 
exist  for  internal  flow  at  temperatures  above  4°C. 

For  data  obtained  in  an  open  channel,  the  first  type  of  correlation  that  one  is  tempted  to  consider 
is  that  of  semi-infinite  flow  over  a  flat  plate.  This  is  an  external  flow  configuration  and  the 
characteristic  length  is  taken  as  thedistance  from  the  leading  edge  of  the  plate.  This  flow  is  developing 
(as  opposed  to  being  fully  developed);  that  is,  the  momentum  or  velocity  and  thermal  boundary  layers 
have  not  reached  theirmaximum  thicknesses.  The  Von  Karrnan  analogy  (Von  Karman  1 939)  is  widely 
used  for  turbulent  flow,  where  the  fluid  properties  are  based  on  the  film  temperatures  and  the  subscript 
x  refers  to  entrance  length: 


2 


(4) 


Nux  = 


0.0292  Pr  Re?'8 


1  +  0.855 {^*r~  1  +  *n  [l  +  5/6  (Pr-  1 )]} 


Gilpin  et  al.  ( 1 980)  conducted  an  experiment  with  ice  on  a  flat  plate  and  developed  a  correlation,  based 
on  boundary  layer  thickness,  which  Lunardini  et  al .  ( 1 986)  converted  to  a  correlation  based  on  entrance 
length,  for  comparison  purposes 

Nux  =  0.00284  Pr  Re?'888  (Pr  =  13).  (5) 

For  laminar  flow,  Chapman  (1967)  suggests 

Nux  =  0.332  Rel72  Pr1/3  .  (6) 

The  velocity  profile  above  a  flat  plate  is  considered  two-dimensional,  while  that  of  the  open  channel 
is  three-dimensional.  The  data  presented  in  this  report  were  obtained  at  the  centerline  of  the  channel, 
and  from  this  point  of  view  it  may  be  acceptable  to  compare  the  data  with  these  correlations.  However, 
a  more  significant  difference  1  ies  in  the  comparison  of  fully  developed  vs  developing  flow,  and  for  this 
reason  the  use  of  the  entrance  length  as  a  parameter  is  given  only  limited  analysis.  Later  in  this  report 
it  is  shown  that  for  these  experiments  the  flow  was  fully  developed. 

Correlations  obtained  from  analysis  of  tube  flow  are  more  appropriate  for  comparison  to  open- 
channel  flow  than  the  flat  plate  correlations  because  of  the  greater  similarity  between  velocity  profiles. 
Vanoni  ( 1941 )  demonstrated  that  the  Prandtl  logarithmic  velocity  distribution  law  for  tube  flow  also 
applies  to  infinitely  wide  open  channels.  For  turbulent,  fully  developed  flow  the  following  equations 
are  often  used.  The  McAdams  (1940)  equation  for  circular  pipe  flow  is 

Nud  =  0.023  Pr0'3  Rep8 .  (7) 

The  subscript  D  refers  to  the  use  of  the  pipe  diameter  as  the  characteristic  length.  For  noncircular  tubes 
the  hydraulic  diameter  (//)  is  often  substituted.  The  equation  is  similar  to  the  Dittus-Boelter  equation 
(Eckert  and  Drake  1959): 

Nu  =  0.0243  Re)}8  Pr0-4 ,  for  heating  (8a) 


Nu  =  0.0265  Rep8  Pr°  3  ,  for  cooling  . 


Petukov  and  Popov  (1963)  suggest 


Nud  = 


RepPr  (Cf/2) _ 

1.07  +  12.7  (Pr2/3  -  1  )  fCrf 2 


where  Cf/2  =  (2.236  In  ReD-4.639)-2. 

Sieder  and  Tate’s  (1936)  correlation  is 

Nud  =  0.027  Rep8  Pr0  33  ( — H — )°' 

'^surface' 


where  p  is  the  dynamic  viscosity,  and  the  material  properties  are  based  on  the  bulk  temperature  rather 
than  the  film  temperature.  For  fully  developed  laminar  flow  in  tubes.  Hausen  ( 1 943)  suggests 


3 


Nud  =  3.66  +  ™**(P!L)  RepPr, 

1  +  0.04  [(D/Z.)  ReDPr] 2/3 


(11) 


where  D  and  L  are  the  pipe  diameter  and  length,  respectively.  The  limiting  Nusselt  value  for  fully 
developed  laminar  flow  in  circular  tubes  is  about  3.66.  The  Reynolds  numbers  for  the  transition  region 
are  generally  taken  as  between  2000  and  10,000  (based  on  hydraulic  diameter). 

Clearly,  the  best  comparisons  of  the  data  would  be  expected  with  correlations  developed  from  open- 
channel  experiments.  However,  there  are  few  of  these  correlations  available.  Ashton  and  Kennedy 
(1972)  correlated  their  data  based  on  water  depth  d  for  water  flowing  over  ice  in  an  open-channel 
flume, 

Nud  =  0.078  Red8  (Pr  -  13).  (12) 

Since  the  width  of  the  flume  was  0.6 1  m,  rewriting  eq  1 2  in  terms  of  hydraulic  diameter  yields 
Nuh  =  0.078  Ren8  (0.61  +  2 cl)-°2  .  (13) 

Lunardini  et  al.  (1986)  correlated  their  open-channel  flow  data  for  water  flowing  over  ice  as 
Nuh  =  0.0078  Pr1/3  Reft927.  (14) 


EXPERIMENTAL  APPARATUS 

In  order  to  obtain  heat  transfer  data  without  ice  that  could  be  compared  with  melting  correlations, 
an  experiment  was  conducted  to  investigate  mechanisms  that  affect  heat  transfer  to  water  flowing  over 
a  chilled  channel  bed.  A  small  experimental  flume  was  designed  and  constructed  in  which  heat  transfer 
coefficients  of  fully  established  flow  could  be  measured  and  observed. 

Flume 

A  15.2-m-long  by  30.48-cm-wide  recirculating  open-channel  flume  was  constructed  in  the  Frost 
Effects  Research  Facility  (FERF)  of  the  U.S.  Army  Cold  Regions  Research  and  Engineering 
Laboratory.  The  flume  was  assembled  primarily  from  plywood  coated  with  fiberglass  epoxy  resin, 

while  the  bed  was  fabricated  from 
ALCOA  type  300  precision  cast  tool 
and  jig  aluminum  plate.  Large  pass, 
Platecoil  heat  transfer  panels  were 
installed  under  the  entire  length  of  the 
aluminum  bed.  Heat  transfer  mastic 
(Tracit  no.  1000)  was  used  to  ensure  a 
continuous  bond  between  the  panels 
and  the  bed.  Figure  1  shows  a  cross 
section  of  the  flume.  The  head  and  tail 
boxes  are  approximately  0.37  m2.  Wire 
mesh  and  ‘'horsehair”  packaging  ma¬ 
terial  were  used  initially  to  reduce 
turbulence  of  the  water  entering  the 
flume  and  a  small  turning  vane  was 
placed  in  the  entrance  to  the  tail  box  in 
order  to  reduce  air  entrainment  in  the 


— *j  30.48cm  |* — 


Figure  I.  Cross  section  of  flume. 
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Tail  Box 


Figure  2.  Schematic  diagram  of  flume. 


water  at  low  flow  depths.  The  sides  and  bottom  of  the  flume  were  insulated  with  Styrofoam  insula¬ 
tion. 

The  design  parameters  of  the  flume  were  developed  from  several  criteria.  A  relatively  small  flume 
was  desired  for  economic  reasons  and  a  smaller  flume  was  expected  to  allow  greater  control  of 
experimental  parameters.  The  desire  to  drill  holes  in  the  edge  of  t  ,e  flume  bed  for  installation  of  ther¬ 
mistors  at  the  centerline  of  the  flume  limited  its  width  because  longer  drill  bits  would  tend  to  drift  ex¬ 
cessively  and  cause  sensor  misplacement  or  break  through  the  surface.  Once  the  width  was  deter¬ 
mined,  other  parameters  were  developed  based  on  the  desire  for  fully  developed  turbulent  flow.  Con¬ 
struction  of  the  flume  and  installation  of  the  instrumentation  took  about  6  months,  with  the  majority 
of  that  time  spent  on  interfacing  and  programing  the  various  instrumentation  components. 

Figure  2  shows  a  schematic  diagram  of  the  flume  and  associated  pumps,  valves  and  flow  meters. 
The  main  pump  was  rated  at  5.6  kW  (7.5  hp)  and  was  capable  of  pumping  3.15x1 0~2m3/s  (500  gpm) 
of  water  through  the  flume.  A  heat  exchanger  located  after  the  pump  discharge  and  in  parallel  with  the 
main  piping  could  be  used  tocontrol  the  temperature  of  the  waterentering  the  flume.  Warm  glycol  was 
pumped  through  this  heat  exchanger  by  pump  no.  1 1  (2.5  x  1 0_3m3/s  at  4 1 3  kPa,  40  gpm  at  60  psi)  of 
the  FERF.  Cold  glycol  was  pumped  through  the  heat  transfer  panels  located  under  the  flume  bed  by 
pump  no.  12(2.5x  10“3m3/s  at  41 3  kPa,  40  gpm  at  60  psi)  and  a  booster  pump  (1.1  kW,  1 .5  hp)  which 
was  installed  in  the  glycol  supply  line. 

Orifice  plates  were  installed  as  indicated  in  Figure  2  and  Tygon  tubing  was  used  to  construct 
manometers,  from  which  flow  rates  were  determined  manually.  The  manometer  calibration  curves  are 
included  in  Richmond  (1988). 
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An  area  approximately  10  m  from  the  entrance  was  designated  tne  test  section.  This  1 .5-m-iong 
section  was  constructed  slightly  differently  from  the  otherparts  of  the  flume,  in  that  here  the  walls  were 
made  of  1 .9-cm-(0.75-in.)  thick  clear  plexiglass.  Additionally,  in  this  section  holes  were  drilled  into 
the  edge  of  the  aluminum  bed  such  that  thermistors  could  be  embedded  along  the  centerline  (Fig.  3). 
By  knowing  the  distance  between  the  thermistors,  the  temperatures  and  the  thermal  conductivity  of 
the  aluminum  plate,  one  could  calculate  the  heat  flow  through  the  bed.  During  assembly,  in  order  to 
obtain  another  measure  of  heat  flow,  heat  flow  sensors  (Micro-Foil,  from  RdF  Corporation,  20453- 
1  SP),  which  had  a  thermocouple  built  into  the  plastic  mount,  were  epoxied  onto  the  underside  of  the 
flume  bed  along  the  centerline,  and  below  the  thermistors  (Fig.  3). 

Thermistors  were  installed  at  each  of  the  locations  indicated  in  Figure  2  in  order  to  obtain  fluid 
temperatures.  For  the  thermistors  installed  within  pipes,  thermistor  wells  were  constructed  from 
0.3 1 75-cm  ( 1  /8-in.)  copper  tubing  by  pinching  and  soldering  one  end  of  the  tubing;  the  tubing  was  then 
bent  so  that  the  tube  would  be  well  into  the  flow.  The  tubing  was  then  installed  within  the  pipe  and  the 
pipe  hole  sealed  around  the  tubing.  All  the  thermistors  installed  in  cavities  were  coated  with  Dow  340 
heat  sink  compound. 

Table  1.  Sensor  nomenclature. 


Data  acquisition 


The  data  acquisition  system  evolved  as  in¬ 

Sensor 

Location 

Accttrac'’1 

formation  was  required  and  associated  instn'  - 

type 

description 

mentation  was  acquired.  Ideally,  all  of  the 
following  data  would  be  collected  at  the  same 
instant: 

Thermistor 

=0.32  cm  below  bed  surface 
at  stations  1-4 

0.02°C 

1 .  Heat  flow  through  the  bed 

2.  Vertical  temperature  profile  in  the  water 

Thermistor 

=  1 .58  cm  below  bed  surface 
at  stations  1-4 

0.02°C 

3.  Vertical  velocity  profile  in  the  water 

4.  Ice  surface  profile  along  the  length  of  the 

Heat  flow3 

Epoxied  to  underside  of  bed 
at  stations  1-4 

test  section  in  the  direction  of  flow  (if  ice  was 
present) 

Thermocouple 

Epoxied  to  underside  of  bed 
at  stations  1-4 

0.1  °C 

5.  Water  depth  and  flow  rates 

Thermocouple4 

Temperature/velocity  probe 

0.04°C 

6.  Water  and  glycol  temperatures 

From  these  data  a  number  of  parameters  were 
to  be  cahulated,  in  particularthe  Nussdt,  Prandtl 

Hot  film  anemometer 

— 

Thermistor 

Water  inlet,  heal  exchanger 

0.02°C 

and  Reynolds  numbers. 

Thermistor 

Water  outlet,  heat  exchanger 

0.02  °C 

Data  acquisition  equipment 

Thermistor 

Airtemperature 

0.02°C 

The  data  acquisition  components  used  in 
these  experiments  are  described  in  Tables  1  and 

2.  Figure  4  shows  a  schematic  of  the  data  acqui¬ 

Thermistor 

Flume  exit,  3  cm  below 
water  surface 

0.02°C 

sition  system.  In  this  instrumentation  assem¬ 

Thermistor 

Glycol  exit  below  flume  bed 

0.02°C 

bly,  an  HP-85  microcomputer  served  as  the 
system  controller,  a  built-in  printer  allowed 

Thermistor 

Glycol  entrance 

0.02°C 

hard  copy  data  output  of  system  temperatures 

Thermistor 

Hume  entrance,  in  head  box 

0.02°C 

and  flow  characteristics  between  tests,  and  a 
larger  external  printer  was  used  to  produce  a 

LVDT6 

Carriage  location 

0.15cm 

hard  copy  of  the  test  data.  Thermistor  resis¬ 
tance  reading .  were  made  and  converted  to  dig- 

LVDT7 

'Estimated 

Ice  depth 

ital  signals  by  the  multiplexer  (as  were  voltage  2Viclory  Engineering  Corp.,  model  T32A1 1 , 2000  ohm 

outputs  from  the  linear  voltage  displacement  3RdF  20453-1  Heat  flow  sensor  w/thermocouple 

transducers  [LVDT],  thermocouples  and  heat  4si  hotnim^nTor^"'3" 

flow  sensors).  The  data  were  then  acquired  by  6CeiescoPT  ioi-60A,60-in.  position  transducer 

the  HP-85.  7Schaevitz  HPD  3000, 3-in.  position  transducer 
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Table  2.  Data  acquisition  hardware. 

Nomenclature  Use 


b.  Underwater  view  of  sensors. 


Figure  5.  Hot  film!  thermocouple  probe. 


An  x-y  carriage  assembly  supported  a  combination  hot-film  sensor  and  28-gauge  copper-constan- 
tan  thermocouple  probe  (Fig.  5).  The  x-y  coordinates  of  the  probe  were  controlled  by  the  HP-85 
computer  via  the  Velmax  motor  controller.  Outputs  from  these  two  sensors  were  input  to  the  HP-85 
via  two  HP-4784  multimeters  (which  could  be  read  by  the  HP-85  faster  than  from  the  multiplexer). 
Velocity  measurements  were  to  have  been  made  with  hot  film  sensors  and  a  TS1  model  1 050  constant 
temperature  anemometer. 
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An  LVDT  could  be  attached  in  place  of  the  hot  film/thermocouple  sensor  in  order  to  profile  the  ice 
thickness  when  present.  If  an  ice  profile  test  were  to  be  run,  the  z  axis  (horizontal  and  parallel  to  the 
flow)  DC  motor  would  be  turned  on  manually  and  the  LVDT  used  to  determine  the  ice  thickness. 

Software 

A  data  acquisition  program  was  written  for  the  HP-85  computer  which  allowed  it  to  communicate 
with  the  peripheral  devices  listed  in  Table  2.  The  program  converted  voltage  and  resistance  measure¬ 
ments  to  appropriate  units,  displayed  or  printed  the  data,  and  then  stored  the  test  data  on  magnetic  data 
cassette  tape.  Figure  6  presents  a  flowchart  of  the  data  acquisition  program.  A  copy  of  the  program  is 
included  as  Appendix  A.  From  Figure  6  it  can  be  seen  that  after  the  program  is  started  a  menu  appears 
which  asks  what  type  of  data  is  to  be  obtained. 

The  five  choices  are  1 )  monitor  any  thermistor  channel,  2)  monitor  system,  3)  conduct  tempera¬ 
ture/velocity  profile  test,  4)  measure  ice  profile,  and  5)  exit  program.  The  selection  of  “monitor  any 
thermistor  channel”  ( 1 )  allows  the  user  to  select  a  desired  thermistor,  the  number  of  readings  to  be 
made  and  the  time  lapse  between  readings.  The  resistance  readings  are  converted  to  temperatures  and 
then  printed  on  the  internal  printer.  “Monitor  system"  (2)  takes  readings  of  the  thermistors,  heat  flow 
sensors,  and  thermocouples,  and  then  does  some  preliminary  analysis.  Typical  output  of  this  routine 
is  shown  in  Figure  7.  The  temperature/velocity  profile  test  routine  (4)  operates  similarly  to  “monitor 
system”  (2)  except  that  at  each  station  measurements  from  the  velocity/temperature  probe  are  made. 
The  program  asks  the  user  for  the  number  of  vertical  positions  (in  the  flow)  and  the  increment  or 
distance  between  readings.  The  number  of  horizontal  (perpendicular  or  transverse  to  the  flow) 
positions  is  entered  along  with  the  distance  from  the  centerline  (+  or  -)  for  each  transverse  position, 
and  the  vertical  locations  are  repeated  at  each  transverse  position.  At  the  end  of  the  measurement 
section  of  this  routine,  basic  data  are  recorded  on  the  magnetic  tape  cassette.  The  data  then  are  printed 
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Figure  6.  Flow  chart  of  data  acquisi¬ 
tion  program. 


on  the  external  printer  along  with  the  results  of  the  initial 
data  analysis  (heat  flow  values,  Reynolds  number,  Nusselt 
number,  etc.).  Appendix  B  contains  a  typical  output  from 
this  routine.  The  ■'monitor”  routines  (1  and  2)  were  used  to 
determine  when  the  system  had  reached  a  steady  state. 

The  ice  profile  routine  allowed  periodic  measurements 
of  an  ice  thickness  on  the  flume  bed  and  allowed  up  to  1 50 
ice  thickness  measurements  to  be  made  over  the  1 ,2-m  test 
section.  Due  to  time  constraints  experiments  with  ice  pres¬ 
ent  were  not  conducted;  hence  this  routine  was  never  used 
and  no  ice  melting  data  were  collected. 

Calibration  and  accuracy 

Tabic  1 ,  as  previously  noted,  contains  the  nomenclature 
for  each  of  the  sensors  used  in  collecting  data.  This  section 
of  the  report  will  describe  how  the  output  of  these  sensors 
was  converted  to  the  desired  units  and  the  associated 
accuracy.  This  section  also  discusses  how  other  values 
needed  to  complete  the  data  analysis  were  obtained. 

Temperature  data  were  obtained  from  two  types  of  sen¬ 
sors,  thermistors  and  thermocouples.  The  thermistors  were 
obtained  with  calibration  data  at  -40°C,  0°C  and  25°C,  the 
resistance  measurements  were  specified  at  ±  0.01  %  and  the 
temperature  at  ±  0.0 1°C.  These  data  were  then  fitted  to  the 
Steinhart-Hart  equation  (Stanley  1967),  which  was  then 
used  to  convert  the  resistance  measurement  to  temperature: 


TIME  11  4 l 25797222 
ft  I R  TEMP  18  98 
PLUME  ENTRANCE  TEMP  4  23 
FLUME  EXIT  TEMP  4  81 
COLD  ClVlOL  ENTR  -33  53 
COLD  GLYCOL  EXIT  -38.83 
HOT  GLYCul  EnTR  14.73 
HOT  GLYCOL  EXIT  14.74 


STATION  l 

HEAT  FLOW  METER  -338.837 
HEAT  FLOW  X  PLATE  <  55H-1686  284 
TEMP  AT  BOTTOM  OF  PLATE-1.64 
THERMISTER  AT  BOTTOM  38 

THERMISTER  AT  TOP  .65 

CALC  PLATE  SURF  .69 


STATION  2 

HEAT  FLOW  METER  -134  234 
HEAT  FLOW  X  PLATE  (..35>-2621  184 
TEMP  AT  BOTTOM  OF  PLATE- . 383 
THERMISTER  AT  BOTTOM  1  89 
THERMISTER  AT  TOP  1  33 
CALC  PLATE  SURF  1  42 


STATION  3 

HEAT  FLOW  METER  -235.175 
HEAT  FLOW  X  PLATE  (. 551-3843  848 
TEMP  AT  BOTTOM  OF  PLhTE-1  491 
THERMISTER  AT  BOTTOM  89 

THERMISTER  AT  TOP  117 
CALC  PLATE  SURF  1  24 


STATION  4 

HEAT  FLOW  METER  -184.335 

heat  flow  x  plate  <  551-3124  333 

TEMP  AT  BOTTOM  OF  PLATE-  478 
THERMISTER  AT  BOTTOM  92 

THERMISTER  AT  TOP  1.83 
CALC  PLATE  SURF  1.28 


Figure  7.  Typical  output  of  “system 
monitor.’’ 


X  =  A  +  B  In  R  +  C  (In/?)3 


(15) 


where  T  =  temperature  in  kelvins 

R  =  resistance  of  the  thermistor 
A,B,C  =  curve  fit  constants 


In  order  to  make  accurate  measurements  from  thermistors  it  is  necessary  to  have  a  low  current  flow 
to  avoid  self-heating.  With  the  HP-342 1 A  this  was  accomplished  by  using  the  300  kilo-ohm  resistance 
range  (10-p.A  current)  during  the  measurements  which  provided  a  measurement  accuracy  of 
±  0.00315%  of  reading  +  3  counts,  in  ohms.  Counts  refer  to  the  value  of  the  last  place  in  the  reading. 
The  accuracy  of  a 2000-ohm  resistance  would  be  3 .06  ohms,  with  the  HP- 342 1 A  at  room  temperature. 
This  corresponds  to  about  ±0.02°C  for  measurements  in  the  0°  to  25°C  range.  At  ambient  temperatures 
different  from  about  23°C,  slightly  higher  inaccuracies  are  possible.  The  thermocouples  used  were 
copper-constantan  and  the  thermocouples  connected  to  the  heat  flow  sensors  had  an  accuracy  of 
±  0.1  °C  when  used  with  the  HP-3421  A.  The  thermocouple  used  to  measure  the  water  temperature 
profiles  was  constructed  from  28-gauge  thermocouple  wires,  carefully  welded  together  and  then 
coated  with  insulating  varnish.  This  thermocouple  was  calibrated  using  a  thermistor  as  a  reference, 
with  each  placed  in  a  constant  temperature  ethanol  bath.  The  reference  junction  for  the  thermocouple 
circuit  was  placed  in  an  ice  bath.  A  new  ice  bath  was  made  at  the  beginning  of  each  day.  A  straight 
line  was  fitted  to  the  resulting  temperature-voltage  data  (Fig.  8);  the  accuracy  of  this  thermocouple 
was  estimated  to  be  0.04°C. 

The  heat  flow  sensors  were  obtained  with  individual  calibration  data  from  the  manufacturer,  and 
no  additional  calibration  was  done.  At  21°C  the  sensors  utilized  were  to  have  an  output  of  6.3442 
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|iV/W-m~2(0.02  |iV/BTU -ft-2-l<r-1 )  and  a  curve  to 
correct  for  different  temperatures  was  also  pro¬ 
vided.  The  output  from  the  multiplexer  was  in 
volts;  therefore,  division  of  the  output  voltage  by 
the  calibration  constant  and  multiplication  by  the 
temperature  factor  should  have  yielded  a  measure 
of  heat  flow.  However,  this  was  not  the  case,  and 
these  sensors  did  not  yield  any  reasonable  data .  The 
only  possible  explanation  for  this  was  that  the 
epoxy  used  to  adhere  these  sensors  to  the  bottom  of 
the  plate  was  acting  as  an  insulator.  This  was 
supported  by  the  thermocouple  output  that  tended 
to  be  lower  than  anticipated  for  the  bottom  surface 
of  the  aluminum  plate,  making  up  the  flume  bed. 

The  constant  temperature  anemometer  was 
obtained  under  a  rental  contract  with  TSI  Inc.  A 
water  probe  calibrator  was  also  obtained  from  TSI  with  the  anticipation  that  frequent  recalibration  of 
the  sensors  would  be  required.  The  flume  was  designed  and  constructed  without  capability  forfiltering 
the  water,  as  we  thought  that  if  the  water  became  too  dirty  it  could  be  drained  and  the  flume  refilled 
with  clean  water.  Cheese  cloth  was  placed  over  a  screen  at  the  entrance  and  used  as  a  coarse  filter. 
(Other  tests  that  were  being  done  in  the  FERF  created  a  dusty  env  i  ronment  and  for  this  reason  the  flume 
was  kept  covered  with  plastic  sheets  which  minimized  dust  settlement  into  the  water.)  It  was 
anticipated  that  dirt  in  the  water  would  have  some  effect  on  the  hot  film  sensors,  and  discussions  with 
a  TSI  representative  indicated  that  the  buildup  of  dirt  would  be  gradual  and  probably  could  be 
accommodated  with  frequent  calibrations.  This  was  not  the  case  and  no  reliable  data  were  obtained. 
The  addition  of  a  diatomaceous  earth  filter  did  not  improve  the  hot  film  sensor  data. 

Since  the  anemometer  did  not  provide  acceptable  flow  velocity  data,  it  was  necessary  to  determine 
the  mean  flow  velocity  from  the  orifice  plates  by  making  manual  measurements  of  the  difference  in 
fluid  height  in  the  manometers.  Of  primary  concern  was  the  orifice  plate  located  in  the  head  box  supply 
line.  Two  different  manometers  were  constructed:  a  manometer  for  high  flow  rates  that  contained 
mercury  and  a  lower  flow  rate  manometer  that  used  water.  Cal  ibration  c  urves  were  produced  for  each 
manometer  and  it  was  estimated  that  the  accuracy  obtained  from  these  measurements  would  be 
±3.15  x  10  nf^/s  (5  gpm).  This  flow  rate  was  converted  to  the  mean  water  velocity  by  using  the 
principle  of  conservation  of  mass. 

The  cross-sectional  area  of  the  flume  occupied  by  the  flow  was  determined  by  reading  the  water 
depth  at  each  test  station  from  a  scale  taped  to  the  plexiglass.  This  depth  measurement  was  accurate 
to  within  0. 1  cm.  The  water  depth  was  multiplied  by  30.48  cm  to  obtain  the  area.  Flow  depths  ranged 
from  10  to  30  cm,  so  the  inaccuracy  of  this  measurement  was  about  ±  1%  in  the  worst  case. 

In  order  to  determine  the  heat  flux  through  the  flume  bed  two  additional  values  were  required,  the 
distance  between  the  thermistors  within  the  plate  and  the  conductivity  of  the  flume  bed  material.  A 
value  of  138  W/m  K  for  the  conductivity  of  the  aluminum  was  obtained  from  the  manufacturer’s 
specification  data  (ALCOA).  An  attempt  was  also  made  to  obtain  the  conductivity  using  CRREL’s 
guarded  hot  plate  apparatus;  however,  no  usable  value  was  obtained  due  to  difficulties  in  using  this 
method  with  high  conductivity  materials.  It  was  therefore  estimated  that  the  conductivity  was  accurate 
to  1 0%  and  that  this  would  include  any  effect  of  the  thermistor  glass  and  the  heat  sink  compound  around 
the  thermistors.  The  holes  were  drilled  into  the  edge  of  the  plate  as  near  to  the  surface  as  possible 
without  risking  breaking  through  the  surface.  The  distance  from  the  surface  to  the  holes  was  later 
checked  with  an  ultrasound  device.  The  average  distance  between  the  thermistors  was  determined  to 
be  0.0127  m  ±0.0001  m;  the  individual  measurements  for  each  test  station  were  used  in  calculating 
the  heat  flux  (Fig.  3). 
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EXPERIMENTAL  PROCEDURE  AND  CONDITIONS 


As  mentioned  earlier  in  this  report,  data  were  collected  only  for  no-ice  conditions.  Data  were 
collected  after  the  system  had  reached  a  steady  state,  identified  by  little  or  nochange  in  the  temperature 
observations  made  using  the  system  monitor  subprogram.  For  high  flow  rate  experiments,  steady  state 
could  be  obtained  in  one  to  two  hours,  while  with  lowerflow  rates  it  could  take  8  hours  or  more  to  obtain 
the  desired  conditions.  Generally,  two  complete  sets  of  data  were  recorded  in  quick  succession  in  order 
to  obtain  eight  data  points  under  nearly  the  same  conditions. 

In  tests  conducted  before  the  water  filter  was  obtained,  the  flume  was  rinsed  out  and  refilled  with 
fresh  water  prior  to  testing  each  day.  After  a  filter  was  obtained,  the  bed  was  wiped  down  and  the  filter 
turned  on.  With  valves  1  and  2  (Fig.  2)  shut  off,  the  main  pump  was  turned  on.  Valves  1  and  2  were 
then  gradually  opened,  and  water  was  added  continually  until  it  was  at  the  desired  depth.  The  water 
depth  was  generally  between  1 0  and  1 5  cm  except  for  some  low  flow  tests  where  the  depth  was  as  high 
as  3 1  cm.  At  this  point  the  manometer  lines  were  purged  of  air  by  manipulation  of  tubing  clamps.  Th~ 
glycol  valves  (4-6)  were  left  open,  pumps  1 1  and  12  were  turned  on,  and  the  temperature  of  the  glycol 
in  these  lines  adjusted.  The  glycol  temperatures  ranged  between  -37°  and-2°C  for  the  cooling  panels. 
Some  experimentation  was  done  in  an  attempt  to  maximize  heat  flow  into  the  glycol  from  the  flume 
bed;  thus  a  flow  rate  of  2.08  x  10'3m3/s  at  -33°C  was  generally  used.  The  warm  glycol  for  the  heat 
exchanger  was  generally  set  at  1 2°C  and  adjusted  as  needed.  The  bulk  water  temperature  of  the  flume 
was  maintained  by  diverting  some  of  the  water  through  the  warm  glycol  heat  exchanger.  The  x-y 
carriage  assembly  was  moved  manually  so  that  the  hot  film  sensor  probe  was  1  mm  (measured  with 
a  “feeler”  gauge)  above  the  plate  surface  centerline.  The  motor  assembly,  when  turned  on,  accepted 
this  position  as  zero  in  the  transverse  (perpendicular  to  flow)  and  vertical  directions. 

The  experimental  parameters  considered  the  most  important  were  the  flow  rate  (Reynolds  number) 
and  bed  surface  temperatures.  Since  most  of  the  tests  were  run  at  the  same  waterdepth  and  similar  bulk 
water  temperatures,  the  Reynolds  number  was  varied  by  adjusting  the  flow  rate.  The  flow  rate  through 
the  flume  was  controlled  by  valves  1  and  2,  and  these  were  adjusted  to  maintain  the  desired  flow  rate 
and  the  temperature  of  the  water.  The  water  temperature  was  adjusted  so  that  the  plate  surface 
temperature  would  be  in  the  desired  range. 

Temperatures  were  monitored  using  the  system  monitor  subroutine  as  described  earlier.  Addition¬ 
ally,  if  the  flume  bed  temperature  dropped  below  0°C,  ice  could  be  observed  on  the  flume  bed  through 
the  plexiglass  walls  of  the  test  section.  The  water  temperature  was  lowered  until  ice  began  to  form  and 
then  raised  slightly.  Adjustments  were  made  as  the  flume  came  toequilibrium.  The  flume  was  assumed 
to  be  at  a  steady-state  condition  when  the  monitored  temperatures  ceased  changing  or  were  changing 
very  slowly  and  randomly. 

Four  flow  conditions  were  of  interest  for  these  experiments:  turbulent  flow  with  the  plate  surface 
temperature  (TJ  near  0°C,  turbulent  flow  with  the  plate  surface  above  4°C.  and  transitional  flow 
(turbulent  to  laminar)  with  Tsnear  0°C  and  above  4°C.  For  the  transitional  or  low  flow  tests,  the  water- 
filled  manometer  was  used  in  place  of  the  mercury-filled  manometer,  allowing  more  accurate  low  flow 
rate  measurements.  Additionally,  lower  flow  rates  (less  than  6.31  x  l0-3m3/s  [100  gpm))  required 
warmer  cold  glycol  temperatures  or  higher  bulk  water  temperatures  in  order  to  prevent  ice  from 
forming  on  the  flume  bed. 


RESULTS  AND  DATA  ANALYSIS 

Data  were  collected  under  the  flow  conditions  described  in  the  previous  section,  and  Table  3 
presents  a  summary  of  the  data.  Entry  numbers  refer  to  the  test  data  obtained  at  each  station,  and  are 
referenced  to  the  data  presented  in  Appendix  C,  which  were  compiled  from  the  data  acquisition  system 
output.  From  the  compiled  data,  Nusselt  number  vs  Reynolds  number  correlations  were  obtained 
based  on  hydraulic  diameter  and  on  entrance  length.  These  calculated  parameters,  along  with  the 
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Table  3.  Data  summary. 


Plate  surface 


calculated  heat  flux,  are  also  listed  in 
Appendix  C,  and  are  cross-referenced 
by  entry  number  to  the  complete  set  of 
data  (Richmond  1 988)  which  are  too  vo¬ 
luminous  to  include  here.  The  data  in 
Appendix  C  were  sorted  into  the  four 
flow  conditions  described  in  Table  3. 

In  order  to  check  that  the  momentum 
boundary  layer  was  indeed  fully  devel¬ 
oped  at  the  test  section,  the  fully  estab¬ 
lished  flow  length  was  calculated  for 
various  flow  rates.  It  was  assumed  that 

the  flow  was  fully  developed  when  the  boundary  thickness  was  the  same  as  the  water  depth.  Using  the 
equation  for  flow  over  a  flat  plate  (Schlichting  1955), 


Flow  condition 

temperature 

Entry  nos. 

Average  Prandll  no.* 

Turbulent 

>4°C 

1-69,232-259 

9.3 

<4°C 

70-183 

10.1 

Transition 

<4°C 

184-215 

5.1 

(turbulent  to 
laminar) 

>4°C 

216-231 

5.3 

*  Based  on  bulk  temperatures. 


8  =  0.37  L 


(16) 


for  the  turbulent  boundary  layer  thickness,  where  is  the  free  stream  velocity,  the  following 
distances  to  attain  established  flow  (()  were  obtained: 


Water  depth  Established* 
Flow  rate  ( GPM )  8  (nt)  flow  length  (m) 


100  0.15  5.81 

300  0.15  7.70 

500  0.15  8.69 


♦Water  temperature  of  10°C. 

The  distance  from  the  inlet  to  the  test  section  was  approximately  1 0  m,  thus  confirming  that  the  flow 
was  fully  developed. 

Figure  9  presents  the  Nusselt  number  vs  Reynolds  number  data  based  on  hydraulic  diameter.  The 
turbulent  data  are  typical  of  those  obtained  for  tubes;  however,  there  is  extensive  scattering  of  the  data 
in  the  transition  region.  The  data  obtained  by  Lunardini  et  al.  (1986),  with  ice,  for  low  Reynolds 
numbers,  fall  in  the  box  shown  in  this  figure. 

The  turbulent  data  for  Ts<  4°C  begins  at  a  Reynolds  number  of  approximately  2  x  104  and  that  for 
T%  >  4°C  at  about  3.5  x  104.  From  this  observation  it  appears  that  the  configuration  with  T$  <  4°C  is 
less  stable  than  that  with  T^>  4°C.  This  is  consistent  with  earlier  observations  suggesting  that  melting 
cases  (Ts  =  Q°)  have  convective  instabilities.  The  low  Reynolds  number  data  seem  to  be  scattered 
around  the  typical  laminar/transition  region  for  nonmelting  systems  rather  than  Lunardini’s  melting 
data,  as  might  have  been  expected  (based  on  temperature  and  density  comparisons).  The  scatter  in  this 
region  was  caused  by  two  problems.  First,  it  was  found  during  calibration  of  the  system  that  cold  spots 
existed  along  the  flume  bed;  these  cold  spots  were  generally  located  near  joints  in  the  heat  transfer 
panels.  In  these  areas  higher  glycol  flow  rates  existed  within  the  panels,  which  subsequently  caused 
higher  heat  fluxes  in  these  same  areas.  Calculations  showed  that  flow  in  the  individual  heat  exchanger 
channels  was  laminar,  thus  heat  transfer  was  limited.  A  small  booster  pump  was  added  to  the  glycol 
supply  line,  but  this  did  not  completely  solve  the  problem.  For  this  reason  different  surface 
temperatures  were  obtained  at  each  station.  This  difference  varied  from  0.2°C  to  1 .0°C.  A  significant 
problem  with  the  data  in  the  low  Reynolds  number  region  is  inaccuracy.  Poor  accuracy  in  the 
calculated  values  (Nusselt,  Reynolds  numberand  heat  flux)  was  caused  by  the  small  differences  in  the 
temperatures  measured  within  the  flume  bed.  An  error  analysis  was  done  following  Holman  (1971), 
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Figure  9.  Nusselt  number  vs  Reynolds  number  (hydraulic  diameter). 


where  the  uncertainty  in  a  function  R  is  defined  by  Wf  and  is  calculated  from  uncertainties  (Wn)  in  the 
independent  variables  Xy  XY  X^ ...,  ^  as  follows: 


Using  this  equation  (App.  D),  the  equations  for  Nusselt  number,  Reynolds  number  and  the  uncertain¬ 
ties  described  above  resulted  in  Figure  1 0,  which  shows  graphically  the  potential  error  for  various  data 
points.  The  largest  source  of  potential  error  was  found  to  be  due  to  small  temperature  differences 
between  the  thermistors  embedded  in  the  flume  bed. 

This  problem  had  been  anticipated  during  design  of  the  flume,  but  was  not  resolved  due  to 
compromises  to  obtain  low  surface  temperatures  at  higher  flow  rates.  If  the  data  had  been  in  the  range 
of  the  low  flow  test  data  of  Lunar- 


dini  et  al.  (1986),  the  potential 
error  in  the  Nusselt  number  would 
have  been  more  tolerable.  This 
can  be  illustrated  by  calculating 
the  maximum  error  associated  with 
the  heat  flux  determination  due  to 
measurement  inaccuracy  in  the 
temperature  alone.  The  error  of 
the  thermistors  was  ±  0.02°C; 
calculating  the  flux  based  on  a 
temperature  difference  of  0.04°C 
(the  maximum  error  associated 
with  the  temperature  difference) 
yields  435  W/m2.  Comparing  this 
value  with  the  average  measured 
value  of  the  low  flow  rate  heat  flux 
(880  W/m2  for  r  >  4°C)  it  can  be 
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Figure  10.  Maximum  potential  error  in  selected  data  points. 
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seen  that  the  maximum  error  is  close  to  50%  at  these  flow  rates.  It  can  be  realized  that  as  the  heat  flux 
increases  the  relative  error  is  reduced.  This  large  error  in  the  low  flow  rate  data  makes  the  quantitative 
values  of  limited  value  and  they  are  not  included  in  any  further  analysis. 

Heat  transfer  correlations  of  the  turbulent  data  in  Figure  9  were  developed  using  a  least-squares 
regression  analysis.  The  equations  were  developed  in  the  form  of 

-MU-  =  A  ReB  (18) 

Pr1/3 

The  use  of  the  Prandtl  number  to  the  1/3  power  comes  from  observation  of  other  correlations  rather 
than  any  clear  trend  in  the  data.  The  correlations  obtained  for  the  turbulent  region  are 

Ts>  4  °C 

NHH  _  0.00128  Reu096  3  x  jo4  <  Re  <  1.5  x  io5  (19) 

Prt/3 

Ts<  4  °C 

=  0.00265  Ren040  1-5  x  104  <  Re  <  1.5  x  105  (20) 

Pr,/3 

where  r.is  the  plate  surface  temperature.  It  can  be  seen  from  these  equations  that  a  higher  (on  the  order 
of  10%)  heat  transfer  rate  is  obtained  for  T,  <  4°C  than  for  T.  >  4°C.  This  is  as  anticipated  from 
observations  of  melting  systems,  which  inherently  include  the  density  inversion.  These  equations  are 
plotted  in  Figure  9  and  10  using  a  Prandtl  number  of  10,  and  bulk  flow  temperature  properties. 

The  data  are  also  shown  in  Figure  1 1 ,  with  parameters  based  on  entrance  length  and  fluid  properties 
based  on  film  temperature.  In  Figure  1 1  it  can  be  seen  that  the  transition  from  laminar  to  turbulent  flow 
occurs  at  about  ReL=  106  as  expected.  It  was  established  earlier  that  the  flow  was  fully  established 
and  the  data  are  presented  in  this  form  only  for  completeness.  The  correlations  obtained  were 


Nux 
Pr  I/3 

=  0.001366  Rex'  071 

Ts  >  4  °C 

(21) 

Nux 

Prl/3 

=  0.00539  Re?'9837 

rs  <  4  °C 

(22) 

Returning  to  the  analysis  based  on 
hydraulic  diameter.  Figure  12  com¬ 
pares  eq  19  and  20  with  some  of  the 
correlations  discussed  earlier.  ^ 

Examination  of  Figure  12  shows  j= 
that  at  low  Reynolds  numbers  there  is,  §  t  o 4 

in  general,  good  agreement  between  5 
all  the  correlations  except  that  of  -£ 

McAdams  (1940,  eq  7).  Both  eq  13  z 
(Ashton  and  Kennedy  1972)  and  eq  1 4  3>  io3 

(Lunardini  et  al.  1986)  compare  very  z 

well  with  eq  9  (Petukov  and  Popov 
1963).  At  higher  Reynolds  numbers, 
the  equations  (eq  1 9  and  20)  obtained  !  0  * 

from  the  test  data  tend  to  diverge  from 
the  other  correlations. 

A  number  of  possible  explanations  Figure  1 1 .  Nusselt  number  vs  Reynolds  number  ( entrance 

for  this  behavior  are  now  examined.  length). 


Reynolds  No.  (entrance  length) 


During  the  experiments,  it  was  thought  that 
this  might  be  caused  by  excess  turbulence 
induced  by  unsettled  entrance  conditions, 
and  because  of  this,  a  different  turbulence- 
reducing  mechanism  was  installed  (hand¬ 
made  aluminum  grid  flow  straighteners). 

The  results  of  additional  tests  (entries  232- 
259)  did  not  significantly  differ  from  ear¬ 
lier  tests.  The  varying  flume  bed  surface 
temperature  could  have  had  an  effect  on  the 
turbulence  in  the  boundary  layer,  but  at 
higher  flow  rates  (i.e.,  higher  Reynolds 
numbers)  this  effect  would  be  reduced. 

Recalling  that  duct  correlations  may  only 
agree  with  wide  shallow  channels  and  also 
noting  that  for  most  of  the  tests  the  width  to 
depth  ratio  was  about  two,  it  is  possible  that 
the  divergence  may  be  caused  by  a  veloc  ity 
profile  effect.  Considering  that  this  effect 
would  be  more  pronounced  at  high  Rey¬ 
nolds  numbers  supports  this  reasoning.  The  other  open-channel  correlations  were  obtained  from 
relatively  wide  shallow  flumes  (width/depth  =  7). 

Further  literature  review  of  the  velocity  distribution  in  narrow  open  channels  vs  wide  channels 
confirmed  that  the  velocity  distribution  in  an  open  channel  is  dependent  on  the  relative  magnitude  of 
channel  width  and  depth  (Rouse  1938).  Since  the  velocity  gradient  at  the  boundaries  will  vary  from 
point  to  point  in  a  cross  section,  then  the  intensity  of  the  boundary  shear  cannot  be  assumed  constant 
over  the  walls  and  bed  of  the  channel.  A  procedure  presented  by  Vanoni  (1977)  provides  a  means  of 
determining  the  average  shear  stress  on  the  bed.  This  is  accomplished  by  separating  the  shear  force 
exerted  on  the  bed  from  that  on  the  walls.  A  bed  shear  velocity  is  then  determined  from  which  a 
corrected  mean  velocity  can  be  obtained. 

The  procedure  followed  was  to  first  determine  the  friction  factor  if)  forthe  flume.  During  initial  tests 
of  the  flume  and  data  acquisition  system,  Manning’s  friction  coefficient  (n)  was  found  to  be  about 
0.012.  The  friction  factor  is  calculated  from 

R-'n  <23) 

where  g  is  the  acceleration  due  to  gravity  and  R  is  the  hydraulic  radius  (area/wetted  perimeter).  The 
value  Re^  was  then  computed.  A  linear  regression  analysis  was  used  to  fit  a  curve  to  data  points 
obtained  from  Vanoni 's  ( 1 977)  plot  which  related  ReH^to/w  (friction  factor  associated  with  the  walls). 
The  equation  developed  was 

u  =  0.253  (ReH//)"°169  .  (24) 

Once /w  was  obtained,  the  friction  factor  due  to  the  bed  ifb)  was  calculated  using 

fb  =  f+  U  (/-/w)  (25) 

b 

where  d  is  the  water  depth  and  b  is  the  bed  width.  Continuing  with  Vanoni ’s  procedure,  the  bed  shear 
velocity  (u»b)  was  then  calculated  using 


Figure  12.  Comparison  of  heat  transfer  correlations 
(Prandtl  number  =  10). 
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u*b  =  ijR^s  = 


(26) 


where  Vavgis  obtained  from  the  mass  flux,  S  is  the  shear  stress  and  /?b  is  the  bed  hydraulic  radius.  The 
mean  velocity  at  the  centerline  was  then  obtained  from  an  equation  presented  by  Yalin  ( 1 977),  which 
relates  the  mean  velocity  to  the  bed  shear  velocity: 

v  =  2.5  «*b  In  (3.32  Re*b)  (27) 

where  Re,bis  the  Reynolds  number  based  on  the  water  depth  and  the  bed  shear  velocity. 

Since  the  procedure  is  intended  to  correct  velocities  in  narrow  flumes,  application  of  the  procedure 
to  wider  flumes,  such  as  that  used  by  Lunardini  et  al .  ( 1 986),  should  have  little  or  no  effect  on  the  mean 
velocity.  The  results  of  calculations  to  check  this  hypothesis  are  shown  below: 

Reynolds  %  change  in  velocity 

number  Narrow  flume  (0.3048  m)*  Wide  flume  (1.2  m)+ 

1.1615  x10s  25%  4.6% 

3.4846  xlO4  29%  8.5% 


*  Width/depth  =  2. 

+  Width/depth  =  7. 

As  can  be  seen,  there  is  only  a  small  effect  on  the  velocity  in  the  wider  flume  at  the  higher  flow  rate, 
while  at  lower  flow  rates  some  effect  is  observed. 

Appendix  C  contains  ReH  values  calculated  using  the  velocity  correction  procedure.  This  velocity 
correction  procedure  is  applicable  only  to  the  turbulent  flow  data,  and  of  course  only  affected  the 
Reynolds  number  calculations.  The  corrected  data  are  plotted  in  Figure  1 3.  Least-squares  regression 
of  this  data  resulted  in  the  following  set  of  equations: 


Nuh  =  0.006  Pr1/3  Ren93 

Ts>  4°C 

(28) 

Nuh  =  0.0142  Pr1/3  Ren864 

Ts  <  4  °C . 

(29) 

The  regression  coefficients  (r2)  were  0.95  and  0.88  for  eq  28  and  29,  respectively.  Figure  14 
presents  the  95%  confidence  intervals  for  the  two  regressions;  these  intervals  only  begin  to  overlap 
at  the  upper  limit  of  the  T$  <  4°C  data,  i.e.,  ReH=  l.lxl  06.  It  may  be  inferred  from  this  analysis  that 
the  two  regressions  represent  the  influence  of  two  different  heat  transfer  phenomena. 

The  only  obvious  difference  between  these  two  data  sets  is  Ts,  which  is  below  4°C  for  eq  29,  and 
which  produces  the  higher  heat  transfer  rate.  In  Figure  1 3  these  equations  are  compared  with  some  of 
the  correlations  presented  earlier.  Equation  29,  which  should  correspond  to  the  melting  correlations, 
falls  slightly  below  those  correlations  (eq  1 3  and  14).  Two  possible  reasons  for  this  difference  are  that 
1 )  other  heat  transfer  mechanisms  are  occurring  in  the  case  of  melting  and  are  enhancing  the  density 
inversion  effect  or  2)  the  average  plate  surface  temperature  was  about  2.0°C  as  opposed  to  0°C  in  the 
melting  case.  Pursuing  this  second  possibility,  we  attempted  to  correlate  the  parameterNuH/(Prl/3Re^J) 
with  the  plate  surface  temperature.  In  Figure  1 5,  the  above  parameter  is  plotted  against  the  difference 
between  4°C  and  the  surface  temperature;  a  least-squares  regression  line  is  fit  to  the  data.  Despite  the 
considerable  scatter  of  the  data  there  is  some  indication  that  the  surface  temperature  (when  below  4°C) 
is  related  to  the  increased  heat  transfer  observed  in  the  correlation.  The  regression  line  nearly  intersects 
the  0.6  value  (from  eq  28)  at  7,= 4  and  the  0.78  value  (from  eq  14,  for  melting  ice)  at  T  =  0  as  might 
be  expected. 
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Nusseit  No./Prondti  No-  (t/3)  Nusselt  No.  (hydraulic  diometer) 


Figure  13.  Nusselt  number  vs  Reynolds 
number  using  the  velocity  correction. 


Figure  14.95%  confidence 
intervals  for  eq  28  and  29. 


In  regard  to  the  comparison  of  eq  29  with  eq  14  (for  melting  ice),  other  flow  characteristics  may  be 
responsible  for  the  increased  heat  transfer  indicated  by  eq  14.  Lunardini  et  al.  (1986)  suggested  that 
othereffects  could  be  the  flux  of  meltwater  injected  at  the  ice/water  interface  and  increased  free  stream 
turbulence  due  to  inherent  waviness  of  a  melting  ice  sheet. 

Gilpin  et  al.  (1978)  conducted  experiments  in  which  the  boundary  layer  temperatures  were 
measured  in  a  developing  flow  of  water  over  heated  and  cooled  plates.  They  observed  similar  convec- 
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Figure  15.  Effect  of  surface  temperature  on  the  correlation  parameter. 


tive  instabilities  in  both  cases  and,  for  the  case  where  the  density  maximum  was  present,  a  roughly 
constant  temperature  (4°C)  for  much  of  the  boundary  layer  thickness  was  observed. 

The  significance  of  the  work  presented  here  is  that  the  density  inversion  and  subsequent  instabilities 
may  in  part  explain  the  increased  heat  transfer  between  melting  and  nonmelting  systems  observed  in 
fully  developed  turbulent  flow. 

It  should  be  noted  that  the  difference  between  eq  28  and  29  is  1 5%  at  a  Reynolds  number  of  5  x  1 04 
and  decreases  at  higher  Reynolds  numbers,  suggesting  a  reduced  effect  of  the  density  inversion  as  the 
Reynolds  number  increases.  This  is  not  supported  by  eq  14  (Lunardini  et  al.  1986);  however,  it  is 
supported  by  eq  13(Ashton  andKennedy  1972)andeq9(Petukov  and  Popov  1963).  It  seems  plausible 
after  examining  the  bed  surface  temperatures  in  Appendix  C  that  this  may  also  result  because  only 
relatively  warm  surface  temperatures  were  obtainable  at  the  higher  Reynolds  number  (thus  skewing 
the  data). 

In  reference  to  Figure  1 3  once  again,  the  difference  between  eq  28  and  7  (McAdams  1 940)  may  be 
attributed  to  inherent  differences  between  open-channel  flow  with  cooling  from  the  bed  as  compared 
with  pipe  flow  with  axisymmetric  cool¬ 
ing.  It  is  interesting  to  note  that  eq  28  has 
nearly  the  same  exponent  as  eq  1 4  (0.927) 
obtained  by  Lunardini  et  al.  (1986). 

Figure  1 6  examines  the  correlations  at 
Prandtl  numbers  other  than  10,  and  it 
should  be  emphasized  here  that  eq  13 
(Ashton  and  Kennedy  1972)  was  based  on 
a  Prandtl  number  of  13  and  is  not  a  vari¬ 
able  in  their  equation.  Trends  similar  to 
those  observed  with  a  Prandtl  number  of 
1 0  are  observed;  however,  this  may  not  be 
significant  due  to  the  limited  Prandtl 
number  data. 

Lunardini  et  al.  ( 1 986)  observed  a  rela¬ 
tionship  between  the  heat  flux  and  Rey¬ 
nolds  number  for  water  flowing  over 
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Figure  16.  Comparison  of  correlations  at  other  Prandtl 
numbers. 
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melting  ice  at  low  flow  rates  (ReH 
<  10,000).  They  observed  a  lim¬ 
iting  heat  flux  of  488.5  W/m2. 
This  phenomenon  is  examined  in 
Figure  17.  Noting  that  only  the 
turbulent  data  uses  the  velocity 
correction,  it  is  seen  that  the  low 
Reynolds  number  data  are  scat¬ 
tered  around  the  value  of  488.5 
W/m2,  albeit  with  more  scatter 
than  shown  by  Lunardini  et  al. 
(1986).  It  is  interesting  to  note 
that  even  with  the  plate  tempera¬ 
ture  above  4°C  (Ts>  4°C)  the  low 
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Figure  17.  Heat  flux  vs  Reynolds  number  (with  velocity 
correction). 


flow  data  scattered  around  this 
same  value. 

Comparison  of  the  Nusselt 
numbersobtainedforthe  low  flow 
conditions  with  T>  4°C  and  7^  <  4°C  indicates  that  there  is  no  effect  of  the  plate  surface  temperature 
on  the  heat  transfer  rate.  This  is  the  opposite  of  what  is  expected  based  on  the  data  of  Lunardini  et  al. 
( 1 986)  and  the  turbulent  flow  data  of  these  experiments.  Lunardini  et  al .  ( 1 986)  reported  an  average 
Nusselt  number  of  75;  in  these  experiments  an  average  value  of  14  was  found  for  flow  in  the  same 
Reynolds  number  region.  Several  differences  exist  between  the  two  experiments;  in  Lunardini  et  al. 
( 1 986)  ice  was  present,  the  bulk  temperature  was  on  the  order  of  8°C  and  the  water  depth  was  about 
0.25  m.  In  the  experiments  reported  here,  no  ice  was  present,  the  bulk  temperature  was  about  30°C  and 
the  water  depth  was  about  0.25  m.  Since  the  lack  of  ice  did  not  significantly  affect  comparisons  for 
the  turbulent  flow  case,  it  is  assumed  that  this  is  not  the  reason  for  the  differences  observed  between 
the  low  flow  rate  data  of  the  two  experiments.  The  remaining  difference  is  the  temperature  profile. 

In  the  case  of  free  convection,  for  enclosed  horizontal  layers,  a  critical  layer  thickness  must  be 
exceeded  before  the  heat  transfer  rate  can  be  enhanced  by  buoyancy  effects.  This  value  is  generally 
described  by  the  Rayleigh  number  and  the  critical  value  is  1 700.  The  Rayleigh  number  is  defined  by 


Ra  =  —  (30) 

Vf 

where  g  is  the  acceleration  due  to  gravity,  P  is  the  thermal  coefficient  of  expansion,  AT  is  the 
temperature  difference,  and  L  is  a  characteristic  length.  Using  this  equation,  with  Ra  set  to  1700  and 
the  upper  and  lower  layer  temperatures  set  at  4°C  and  0°C  respectively,  yields  a  critical  thickness  of 
about  7  mm.  Layer  thicknesses  less  then  7  mm  will  act  as  conductive  layers  (Nu  =  1 ),  while  for  thicker 
layers,  Bernard  cells  will  develop  ;.nd  enhance  heat  transfer. 

A  possible  reason  for  the  increased  heat  transfer  observed  by  Lunardini  et  al.  (1986)  couH  be  the 
larger  distance  between  the  0°  and  4°C  isotherms  compared  to  the  data  of  this  report.  A  review  of  the 
temperature  profiles  for  these  data  in  Richmond  (1988)  showed  that  in  all  cases  the  4°C  isotherm  was 
less  than  1  mm  thick.  The  temperature  profile  data  of  Lunardini  et  al.  ( 1 986)  have  not  been  published. 
Gilpin  et  al.  ( 1978)  examined  the  external  flow  case  of  flow  developing  on  a  flat  horizontal  plate,  with 
the  density  inversion  of  water  present,  and  were  able  to  visualize  an  unstable  region  of  vortex  rolls 
between  the  0°  and  4°C  isotherms.  These  ideas  seem  to  explain  the  apparent  differences  between  the 
two  sets  of  open-channel  data. 
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SUMMARY 


A  small  cross-sectional  area  flume  was  constructed  in  which  both  the  bed  and  entrance  temperatures 
could  be  controlled.  The  flume  was  designed  so  that  thermistors  installed  in  the  bed  allowed  the  heat 
flux  across  'he  bed  to  be  measured.  Flow  rates  up  to  3.15  x  10_2m3/s  (500  gpm)  were  possible.  Row 
velocities  were  determined  from  manometers  measuring  flow  in  the  supply  pipes,  and  a  velocity 
correction  procedure  was  used  to  compensate  forthe  narrow  flume.  Experiments  were  conducted  with 
Reynolds  numbers  based  on  hydraulic  diameter  up  to  2.3  x  105.  The  flume  had  several  limitations.  It 
was  not  possible  to  measure  small  heat  fluxes  across  the  flume  bed  due  to  difficulties  in  obtaining 
temperatures  with  accuracies  greater  than  0.02°C,  and  uneven,  limited  cooling  of  the  flume  bed 
reduced  its  usefulness  at  high  flow  rates.  These  limitations  mainly  affected  the  accuracy  of  the  Nusselt 
and  Reynolds  numbers  at  low  flow  rates  as  indicated  in  Figure  10,  thus  limiting  the  quantitative 
usefulness  of  the  low  flow  rate  data. 

A  slightly  different  heat  transfer  correlation  was  obtained  for  the  flume  as  compared  to  that 
associated  with  pipe  flow,  which  must  be  due  to  differences  in  the  boundary  conditions  between  the 
chilled-bed  flume  and  axisymmetric  pipe  cooling. 

Comparison  ot  the  two  heat  transfer  correlations  (eq  28  and  29)  obtained  from  these  experiments 
showed  that  an  increase  in  heat  transfer  occurs  when  the  plate  surface  temperature  is  below  4°C,  which 
could  be  due  to  the  density  inversion  of  water  at  4°C  and  subsequent  flow  instability.  The  heat  transfer 
rate  appears  to  be  a  function  of  4°C  minus  the  plate  surface  temperature,  which  supports  the  above 
conclusion  that,  as  the  plate  temperature  decreases  below  4°C,  the  heat  transfer  is  nhanced  by  the 
increasing  buoyancy  effect  of  the  density  inversion. 

Extending  these  findings  to  comparisons  of  the  correlations  obtained  by  others  from  water  flowing 
over  melting  ice  indicates  that  the  density  inversion  of  water  plays  an  important  role  in  the  increased 
heat  transfer  rate  in  fully  ueveloped  turbulent  flow,  although  other  effects  are  not  ruled  out.  Other 
researchers,  (e.g.,  Gilpin  et  al.  1978)  have  shown  that  cooling  is  similar  to  heating  over  a  flat  plate  if 
the  density  maximum  exists;  that  is,  if  the  plate  surface  temperature  is  below  4°C  in  the  entrance  region. 
This  work  demonstrates  that  the  instabilities  continue  into  the  fully  developed  flow  regime  and 
explains  in  part  the  increased  heat  transfer  observed  in  fully  developed  melting  systems. 
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APPENDIX  A:  DATA  ACQUISITION  PROGRAM 


10  DISP  "FLUME  DATA  COLLECTION" 

Z0  ! 

30  CLEAR  7 
40  DISP 

50  OISP  "SET  TIME  BY  TYPING  SETTIME  (HOUR  OF  DAY  *3600  +MIN  OF  HOUR  *6 
0 , ( DAY  OF  MONTH)" 

60  DISP  "EX.  SETTIME  8*3600+5*60,12" 

70  DISP  "PRESS  CONT" 

80  PAUSE 
90  I 

200  OPTION  BASE  1 
210  IMAGE  4D.DD 
220  IMAGE  5D.DDD 

230  DIM  R(  8 ) ,P(40) , S(  20, 3) , A( 7 ) ,P1  (  50) ,PZ(  50) ,P3( 50) ,P4<  50) ,PS( 50) 

235  DIM  L 1 <  4) ,LZ( 4  > , T6<  4),N1(4),NZ(4) ,N4( 4) ,N5<  4) ,R5( 4) ,R6( 4) . T( 4 ) 

240  MAT  READ  S 

250  I  * . * . *********** 

260  DEF  FNV(T)  -  ( -(  . 000004 37S*T* 3 >  +  . 001 3*T*Z-. 14B8*T+7. 9 )».0004 1 339  ! 
VISCOSITY  OF  WATER 

270  DEF  FNK(T)  -  (-(  .000001 1 1 *TAZ )+. 0006889 *T+. 306)* 1 .7303  I  COND.  OF 
UATER 

280  DEF  FNP(T)  -  -< .00001563*TA3)+.0046ZS*TAZ-.S178*T+Z5.6  1  Pr  no 
290  DEF  FND(T)  -  <6Z.23613+.01Z0Z7*T-.000ZZ038*TAZ+.0000007951*TA3)*16 
.018  1  DEN  OF  UATER>4 

300  DEF  FNDKT)  -  (  .99987+. 000006749996*T-. 00000541 659*TA2-. 0000025000 
35*  T A  3+4 . 1 667 1 67E-7*  T  A  4 ) *  1 000 

310  OEF  FNDZ(T)  -  1004. 91Z359*( 1 .088-. 001*T)  !  DEN  OF  GYLCOL 
320  DEF  FNCZ(T)  -  1000*( Z.961+.005*T)  I  Cp  OF  GYLCOL 

330  I  . . ******* 

340  DISP  "UHAT  TYPE  OF  TEST  IS  THIS?" 

350  DISP  "TEMP/VEL  PROFILE- 1  " 

360  DISP  "ICE  PROFILE-Z" 

390  DISP  "MONITOR  SYSTEM  -5" 

400  DISP  "MONITOR  ANY  THERMISTER  CHANNEL  -6" 

405  DISP  "EXIT-7" 

410  INPUT  Z 
420  PRINTER  IS  2 
430  I 
435  K-0 

450  IF  Z-5  THEN  750 
460  IF  Z-6  THEN  4180 
465  IF  Z-7  THEN  STOP 

468  OUTPUT  707 

469  PRINTER  IS  707.80 

470  !  *•**•**•*• . ••••••*• 

471  I  LOAD  8300  PROGRAM 

472  GOSUB  6700 

473  UAIT  1000 

474  CONTROL  10,9  i  137 
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475  WAIT  1000 

480  I  ENTER  TEST  IDENTIFICATION  DATA 
490  I 

500  DISP  "ENTER  TEST  NO.  ( _ >"  0  INPUT  T# 

510  CREATE  TS.Z0 
520  ASSIGN#  1  TO  T* 

530  DISP  "ENTER  DATE  (_  MON  85)" 

540  INPUT  D$ 

550  DISP  "ENTER  UATER  DEPTHS  .cm" 

560  DISP  "STATION  1  "  6  INPUT  D1 < 1 > 

570  DISP  "STATION  2  "  0  INPUT  DKZ) 

580  DISP  "STATION  3  "  0  INPUT  DK3) 

590  DISP  "STATION  4  “  0  INPUT  DK4> 

600  I  ENTRANCE  LENGTHS  .METERS 

610  L5( 1 >“10.55 

620  L5( 2 )“10. 86 

630  L5<  3)“1 1 . 164 

640  LS<4)-1 1.465 

650  DISP  "ENTER  FLOU  RATES  .gal/min" 

660  DISP  “UATER"  0  INPUT  FI 
670  DISP  "COLD  GLYCOL"  0  INPUT  F2 
680  DISP  "HOT  GYLCOL"  0  INPUT  F3 
690  DISP  “HOT  WATER"  0  INPUT  F4 
700  I 

710  !  SETS  UP  CALIBRATION  CONSTANTS  FOR  HOT  FILM 

711  I  SENSOR  #05 

71 Z  DISP  "DATA  IN  PROGRAM  IS  FOR  SENSOR  #05" 

713  DISP  "INPUT  RESISTANCE  FOUND  FROM  UATER  TEMP"  0  INPUT  T8 

714  T8— 546.4S3+T8*98.00Z  I  Tef  DEG  C 

715  DISP  "INPUT  THE  RESISTANCE  SETTING"  0  INPUT  R90  R8-R9+. Z6+40 

716  T9— S46.453+R9*98.00Z  I  Tsf  DEG  C 

717  F9-T9* 1 . 8+3Z 

718  P9-6.48/R9*(  R8/46.74),‘Z*(  (  T9-T81/61 ,74)*( FNK( F9  )/FNK(  191.481) 

719  F7-<T9+T8)/Z«1.8+3Z 

7Z0  Q9“P9»(FNV( 135.91 )/FND< 135.91 >/( FNV( F7)/FND( F7) ) ) * . 4* 152 .09033 

721  P9-P9*36.4665 

722  DISP  "P9-"iP9."09-"j09 
7Z3  DISP  "VISCOSITY" »FNV(F7) 

730  I  . . . . . 

740  ! 

750  I  NO  ICE  TEST  OR  MONITOR 
760  I  SYSTEM 

770  I 

840  Ll< 1 )“Z.751Z1 
850  LI (Z 1-4.7049 
860  LK31-6.657Z 
870  LI ( 4)-8. 61 Z9 

875  LZ< 1 )«.01Z8Z7 

876  LZ< Z)-.01ZGZ 

877  LZ( 3)-. 01 Z7 
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878  L2<4>-. 01267 
880  I 

890  I  READ  SYSTEM  DATA  CHANNELS 
900  ! 

910  I 

920  T-TIME/3600  <  TIME  IN  HRS 
930  OUTPUT  709  !,•0PN•, 

940  OUTPUT  709  s"CLS00" 

960  OUTPUT  709  ; "F3N5RA0R5LS2-8; T3" 

970  ENTER  709  »  A(  1 ) . A(  2 ) , A(  3) , A(  4 ) , A( S) ,A( 6> ,A(  7) 

980  OUTPUT  709  •,,,OPN" 

990  IF  Z-Z  THEN  2120 

1000  I 

1001  FOR  1-1  TO  4 

1002  H-19+I 

1003  C-23+I 

1004  R1-2*I 

1005  R2-R1+1 

1010  I  TAKE  DATA  AT  EACH  STATION 
1020  I 
1030  t 

1040  OUTPUT  709  i "OPN" 

1050  OUTPUT  709  ; "DCVM i H 
1060  ENTER  709  t  H< I ) 

1070  OUTPUT  709  { "TEM“ ; C 
1080  ENTER  709  i  T(  I ) 

1090  OUTPUT  709  i "OPN" 

1100  IF  I>1  THEN  1170 
1110  OUTPUT  709  5"OPN" 

1120  OUTPUT  709  ?“CLS01" 

1130  OUTPUT  709  i "LSZ-9" 

1140  OUTPUT  709  ; "F3N5RA0R5T3" 

1150  ENTER  709  j  R(  1  ) , R(  2  )  ,R(  3)  ,R< 4 >  ,R(  5>  ,R(  6)  ,R( 7)  ,R(  8) 

1160  OUTPUT  709  i "OPN” 

1170  IF  Z<>5  THEN  GOSUB  Z560  I  SUBROUTINE  TO  MOVE  CARRIAGE 
1180  I 
1190  I 
1200  ! 

1210  ! 

1220  I 

1230  !  PROBE  TEMPERATURES 

1240  IF  Z-S  THEN  GOSUB  2750  fi  GOTO  1640  I  SKIPS  PROBE  PROFILE 
1260  IF  I>1  THEN  1510 
1430  CONTROL  10,9  i  136 

1440  OISP  "NUMBER  OF  POINTS  <Y  POSITIONS)" 

1450  INPUT  13 

1460  DISP  "INCREMENT.  +  IS  UPWARD" 

1470  INPUT  14 

1480  DISP  "#  OF  X-POSITIONS" 

1490  INPUT  IS 
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1500  REDIM  P1(I3),PZ(I3),P3(I3) ,P4< 13) ,P5( 13 > 

1510  FOR  K-1  TO  15 
1520  DISP  "X-POSITION" 

1530  INPUT  X 

1550  OUTPUT  10  “RUN  200" 

1560  UAIT  500 

1570  CONTROL  10,9  i  141 

1580  UAIT  1000 

1590  OUTPUT  10  USING  "K ,K ,K ,K ,K"  i  Xj  "  . '' ;  I4»  -  , "  5 13 

1600  UAIT  1000 

1610  CONTROL  10,9  t  141 

1620  60SUB  1750 

1630  NEXT  K 

1640  NEXT  I 

1650  IF  Z-1  THEN  60SUB  6210 
1655  IF  Z-1  THEN  ASSIGN#  1  TO  * 

1660  GOTO  340 

1740  I  .##•••####•######**##*## 

1745  I  SUBROUTINE  FOR  PROBE  MEASUREMENTS 
1750  FOR  J-1  TO  13 
1760  UAIT  6000 
1770  P1(J)-I 

1775  DISP  "NOU  READING  X-Y  COORD” 

1780  ENTER  10  USING  "#,K“  1  P2<J>#  ENTER  10  USING  "#.K"  <  P3(J> 

1785  P3<  J>— <  1*P3<  J) )+.  1 

1786  P5<  J)-0 

1790  FOR  K9-1  TO  10 

1800  OUTPUT  703  j "F 1 R-2RAZ 1 NST3” 

1820  ENTER  703  i  PG  I  VELOCITY 

1821  P5( J)-P5( J)+P6 

1822  NEXT  K9 

1823  P5( J)-P5( J)/ 10 

1827  OUTPUT  723  |”F1R-2RAZ1NST3" 

1828  ENTER  723  »  P4<J) 

1840  CONTROL  10.9  1  141 

1850  OUTPUT  10  rCONT"  6  UAIT  1000 
1860  NEXT  J 
1870  GOSUB  2750 
1880  RETURN 

2110  l  . . . . . 

2120  I  SUBPROGRAM  TO  PROFILE  ICE 
2130  DISP  "ICE  PROFILE  SUBROUTINE" 

2140  DISP  "INPUT  NUMBER  OF  READINGS" 

2150  DISP  "DETAILED  -  150" 

2160  DISP  "6R0SS  -  30" 

2170  DISP  "TEST  -  5"  8  INPUT  J9 

2230  OUTPUT  709  i"CLS01" 

2240  OUTPUT  709  »"LS2-9" 

2250  OUTPUT  709  1 "F3N5RA0R5T3" 

2270  ENTER  709  i  R<  1 )  ,R<  2 >  ,R(  3)  ,R(4>  ,R<5>  ,R<  6)  ,R(7)  ,R< 8) 
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2275  OUTPUT  709  j"OPN" 

2280  OUTPUT  709  j  ,'OCV20-23“ 

2290  ENTER  709  »  H( 1) ,H( 2 > ,H< 3> ,H(4 > 

2300  OUTPUT  709  j”TEM24-27" 

2320  ENTER  709  j  T( 1 ) ,T< 2 > ,T( 3) ,T< 4) 

2325  OUTPUT  709  i"OPN“ 

2330  BEEP 

2340  01 SP  " CARRIAGE  SHOULD  BE  AT  POINT  1  NOW 

2345  01 SP  M TURN  ON  CARRIAGE  NOW 

2350  01 SP  "PRESS  CONT  WHEN  READY"  6  PAUSE 

2360  FOR  1-1  TO  J9 

2370  OUTPUT  709  r,DCV28,29" 

2380  ENTER  709  j  P1<I),P2<I)  !  PI -CARRAGE .P2-ICE  DEPTH 
2390  BEEP 
2400  WAIT  100 
2405  NEXT  I 

2410  01 SP  "PAUSED"  8  PAUSE 
2420  I  CONVERT  TEMP  DATA 
2430  FOR  1-1  TO  4 
2440  GOSUB  2750 
2450  NEXT  I 

2460  I  CONVERT  PROFILE  DATA 

2470  FOR  1-1  TO  J9 

2480  PI (I )-P1( I >*15.54467  !  c« 

2490  P2( I >-P2( I >*( 10/(3*2.54) >  I  cm 
2500  PRINT  I ,P1 ( I > ,P2( I ) 

2510  NEXT  I 
2520  GOSUB  6220 
2530  GOTO  340 

2550  !  . . . . . 

2560  !  SUB  PROGRAM  TO  MOVE  CARRIAGE 

2590  DISP  "MOVE  CARRIAGE  TO  STATION  ";I  8  BEEP  8  WAIT  100  8  BEEP 
2600  DISP  "ONE  BEEP  AND  CARRAIGE  IS  IN  CORRECT  PLACE" 

2610  DISP  LKI) 

2620  WAIT  8000 

2630  OUTPUT  709  ?"DCV28" 

2640  ENTER  709  i  L 
2650  DISP  L 

2660  IF  L<L1 ( I >+.05  AND  L>L1(I>-.05  THEN  2680  ELSE  2670 

2670  BEEP  8  WAIT  100  8  BEEP  8  WAIT  100  8  BEEP  8  GOTO  2620 

2680  BEEP  150,300  8  DISP  "AT  POSITION  "»I  8  OUTPUT  709  »"OPN"  8  RETURN 

2720  !  *********************** 

2750  I  SUBPROGRAM  TO  CONVERT  DATA 
2760  DISP  "CONVERTING  DATA" 

2770  !  CONVERTS  DATA  AND  WRITES  TO  BUFFER 
2775  IF  K>1  THEN  3800 
2780  IF  I>1  THEN  3050 
2790  FOR  J1-1  TO  5 

2800  A< J1 >-1/(S( J1 .1 )+S< J1 ,2>*L0G(A< J1 >>+S< J1 . 3>*L0G< At J1 > >* 3>-273. 15 
2820  NEXT  J1 
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2830  A(G>-1/<S< 1S,1 >+S( 16, 2 )*LOG( A<  6) >+S( 16,3>*L0G< A(G> >*3>-273. 15 
2840  A( 7>-1 /<  S( 17,1 )+S(  17,2 )*LOG< A< 7 > >+S(  17,3)*L0G( A(7) )*3)-273. 15 
2860  PRINT  "TIME" ; T 
2870  PRINT  “AIR  TEMP"s 
2880  PRINT  USING  210  ;  A<5> 

2890  PRINT  "FLUME  ENTRANCE  TEW5" ; 

2900  PRINT  USING  210  ;  A<  1  > 

2910  PRINT  "FLUME  EXIT  TEMP"; 

2920  PRINT  USING  210  ;  A< 2 ) 

2930  PRINT  "COLD  GLYCOL  ENTR"; 

2940  PRINT  USING  210  ;  A<4) 

2950  PRINT  "COLD  GLYCOL  EXIT"; 

2960  PRINT  USING  210  ;  A<3> 

2970  PRINT  "HOT  GLYCOL  ENTR"; 

2980  PRINT  USIN6  210  ;  A(6> 

2990  PRINT  "HOT  6LYC0L  EXIT"; 

3000  PRINT  USING  210  ;  A(7> 

3010  PRINT 

3020  IF  205  THEN  PRINT#  1  ;  A<1  > , A<  2 ) , A<  3> , A<4> , A<  5) . A<  G> . A< 7> ,  13 
3030  I 

3040  11-7  6  12-8 

3050  R(  R 1  —  1  )-1/<S(I1 ,1 )+S(I1 ,2)*L0G(R(R1-1 )>+S(I1 , 3)*L0G(  R(  R1 - 1 >)“3)-273. 15 
3060  R(  R2- 1 )-1/(S(I2,1 >+S(  12,2 >*L06( R(  R2- 1 ) )+S(  12 , 3)*L0G( R( R2- 1 ) ) A3>-273. 1 5 
3070  H( I )-H( I >*50000000* (-( .0014*(T(  I >*1.8+32 >>  +  1.1 > 

3080  H< I )-H( I )/3. 152  !  <U/nA2> 

3090  K 1-1 38  !  U/nC 

3100  HI  ( I )-K 1 *( R( R2- 1 >-R(  R1 - 1 ) )/L2< I >  !  U/n"2 

3210  T6< I )«< R< R1- 1 )-R< R2- 1 ) )*(  1 .9144- 1 . 5893 )/(  1 . 589 3- . 3 1 864 >  +R( R I -  1 ) 

3220  IF  Z-5  THEN  GOTO  3510 
3230  IF  I>1  THEN  GOTO  3510 
3240  I  ENERGY  ANALYSIS 
3250  I  HEAD  BOX  TO  TAIL  BOX 
3260  T1 -A( 2 )-A< 1 > 

3270  IF  A(1>>4  THEN  PI -FND< ( A< 1 >+A(  Z > >/Z* 1 . 8+32 > 

3280  IF  A( 1 ><-4  THEN  P1-FND1 ( ( A( 1 >+A( 2 > >/Z > 

3290  SI -4 186  !  Cp  OF  WATER 

3300  FI -FI *.0000630902  I  FLOW  RATE  IN  k»'3/SEC 
3310  E1-S1*P1*T1*F1*.001  I  KWATTS 
3320  I  COLD  GLYCOL  ENERGY 
3330  TZ-A( 3)-A(4) 

3340  SZ-FNCZ((A(4)+A(3)>/2> 

3350  PZ-FNDZ( ( A<  4 >+A<  3) >/2 ) 

3360  F2-F2* .00006 3090Z 
3370  E2-SZ*P2*TZ*F2* .001 
3380  !  HOT  GLYCOL  ENERGY 
3390  T3-A( 6>-A(7> 

3400  S3-FNCZ<  <  A<  6)+A(  7 > )/2  > 

3410  P3-FNDZ( ( A<  6>+A<  7>  >/2  > 

3420  F3-F3*. 0000630902 
3430  E3-S3*P3*T3*F3*  .001 
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344®  IMAGE  40.000 
345®  PRINT  "FLUME  ENERGY 
346®  PRINT  USING  344®  ;  El 
3470  PRINT  "COLD  GYLCOL  EN" s 
348®  PRINT  USING  3440  i  EZ 
3490  PRINT  "HOT  GYLCOL  EN"; 

350®  PRINT  USING  3440  ;  E3 

3510  IF  Z<>5  THEN  PRINT#  1  »  R( R1 -  1 ) ,R( RZ- 1 ) ,H< I > , T< I > 

35Z0  I1-I1+Z  C  IZ-IZ+Z 
353®  I 
354®  I 

3550  PRINT 

3551  IF  Z-1  AND  1-1  THEN  OUTPUT  707 
355Z  IF  Z-1  AND  1-3  THEN  OUTPUT  707 
3555  IF  Z-1  THEN  GOTO  376® 

356®  PRINT  "STATION  ";I 
3570  PRINT  "HEAT  FLOW  METER" j 
3580  PRINT  USING  ZZ0  j  H(  I ) 

3590  PRINT  "HEAT  FLOW  X  PLATE  (.55)"; 

3600  PRINT  USING  ZZ0  j  H1(I> 

3630  PRINT  “TEMP  AT  BOTTOM  OF  PLATE"!T<I> 

3640  PRINT  "THERMISTER  AT  BOTTOM" i 
3650  PRINT  USING  Z1®  i  R(RZ-I) 

3660  PRINT  "THERMISTER  AT  TOP  "i 
3670  PRINT  USING  Z10  j  R<R1-1) 

368®  IMAGE  30. 3D 

3730  PRINT  "CALC  PLATE  SURF"j 

3740  PRINT  USING  Z10  j  T6(  I ) 

3750  PRINT 

3760  IF  Z-5  THEN  4140 
3770  IF  Z-Z  THEN  4140 
3780  PRINT 
3785  PRINT#  1  i  13 
3800  FOR  JZ-1  TO  13 
3805  PRINT#  1  j  P5(JZ) 

381®  P4(  JZ )-. 1 3573+Z5. 44371 *( P4<  JZ )» 1000) 

3815  P5( JZ )-( ( P5( JZ  >* Z-P9  >/Q9 )*Z . 5* . 3048 

38Z0  NEXT  JZ 

38Z 1  PRINTER  IS  707.80 

38ZZ  OUTPUT  707  i"" 

38Z 3  I MAGE  SX . "STATI ON" . 6X , " XCOORD" . 5X . " YCOORD" . 4X , " TEMP( C ) “ . 4X . " VEL (  M/S  > " 
38Z4  PRINT  USING  38Z3 

38Z9  IMAGE  7X.D, 10X.ZD. ZD.6X, ZD. ZD.5X ,ZD. ZD.7X.ZD. ZD 

3830  FOR  M-1  TO  13 

3831  IMAGE  7X.D, 10XfZD.ZD.6X,ZO.ZD,5X,ZD.ZD,7X,ZD.ZD 
383Z  PRINT  USING  38Z9  j  PI  (  M) ,PZ<  M) ,P3( M) ,P4(M) ,P5( M> 

3835  PRINT#  1  »  PI ( M> ,PZ( M) ,P3( M) ,P4( M> ,P5( M> 

3837  NEXT  M 

3838  PRINT  6  PRINT  C  PRINT 

3839  IF  K>1  THEN  4140 


29 


3840  I  CALCULATE  Re  AND  Nu 
3850  » 

3850  81 -P4(  3 )* 1 . 8+3Z  I  BULK  TEMP 
3870  BZ-T6< I >*1 .8+32  !  SURFACE  TEMP 
3880  B3-<  BJ  +BZ  )/Z  I  FILM  TEMP 
3890  B9*(  B3-3Z  )/ 1 . 8 

3900  IF  B9<-4  THEN  B4-FNDKB9)  I  DEN  WATER  <-  4 
3910  IF  B9>4  THEN  B4-FND(B3)  I  DEN  WATER  >4 
3920  B8-FNV(B3)/B4  !  DYNAMIC  VI SCSI TY 
3930  B5-FNK(B3) 

3940  B7«FNP(B3>  !  PRANDTL  NO 
3950  H 1 *H< I ) / ( T6( I > -P4( 3 ) ) 

3960  HZ-H1<I)/<T6(I>-P4(3>) 

3980  Vl( I )-FJ/( 30.48*D1( I )/( 100*100)) 

3990  H9-4* 30.48*01 <I)/(Z*D1(I>+30.48>/100 

4000  R5<I)-V1(I)«H9/B8  i  Re.Hd 

4010  R6<I)-V1<I)*L5(I)/B8  I  Re,L 

4020  Ni< I )-H1*H9/BS  I  Nu, Hd 

4030  NZ( I )-HZ*H9/BS 

4050  N4(I)-H1*L5(I)/BS 

4060  N5(I)-HZ*L5(I)/BS 

4080  PRINT  "K  OF  WATER" ;B5 

4090  PRINT  "HEAT  TRAN  CO"iHZ 

4100  PRINT  "HYDR  DIAM";H9 

4110  PRINT  "VISCOSITY  OF  WATER" ; B8 

41 Z0  PRINT  "VELOCITY  OF  UATER"5V1<I> 

4130  PRINT  "DENSITY  OF  WATER  "|B4 

4140  RETURN 

4150  I 

4160  I 

4170  ! 

4180  I  MONITORS  ANY  THERMISTER  CHANNEL 
4190  CLEAR  709 

4200  DISP  "WHICH  CHANNEL?"  0  INPUT  G 
4210  DISP  "NUMBER  OF  READINGS?"  •  INPUT  G1 

4220  DISP  "TIME  BETWEEN  READINGS  (SEC)?"  «  INPUT  GZB  GZ»GZ*1000 
4230  IF  GO)  THEN  GOTO  4340 
4240  FOR  1-1  TO  G1 
4250  OUTPUT  709  j"CLS0" 

4260  OUTPUT  709  < "F3N5RA0R5LSZ i T3" 

4270  ENTER  709  j  G3 

4280  G4-1 /(  S(  1 ,1 )+S< 1 ,2 )*LOG( G3)+S( 1 , 3>*L0G(G3)*3>-273. 15 
4290  OUTPUT  709  i"OPN" 

4300  PRINT  USING  210  <  G4 
4310  WAIT  GZ 
4320  NEXT  I 
4330  GOTO  350 

4340  IF  GOZ  THEN  GOTO  4450 
4350  FOR  1-1  TO  61 
4360  OUTPUT  709  |"CLS0" 
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4370  OUTPUT  709  i "F3NSRflOR5LS3i T3“ 

4380  ENTER  709  i  G3 
4390  OUTPUT  709  rOPN" 

4400  G4-1 /( S(Z,1 >+S(Z,Z )*LOG( 63 >+S( 2 , 3>*LQG( G3)A  3)-Z73. 15 

4410  PRINT  USING  Z10  s  G4 

4420  WRIT  GZ 

4430  NEXT  I 

4440  GOTO  350 

4450  IF  G<>3  THEN  GOTO  4550 
4460  FOR  1-1  TO  G1 
4470  OUTPUT  709  ;"CLS0" 

4480  OUTPUT  709  ; "F3N5RR0R5LS4; T3“ 

4490  ENTER  709  i  G3 
4500  OUTPUT  901  !"0PNM 

4510  G4-1/<  S( 3, 1 )+S( 3,Z  >*L0G<G3)+S( 3,3)*L0G(G3)*3)-Z73. 15 
45Z0  PRINT  USING  Z10  ;  G4 
4530  URIT  GZ  6  NEXT  I 
4540  GOTO  350 

4550  IF  G<>4  THEN  GOTO  4GG0 
4560  FOR  1-1  TO  G1 
4570  OUTPUT  709  j"CLSO" 

4580  OUTPUT  709  i,‘LS5" 

4590  OUTPUT  709  | "F3NSRR0R5LS5? T3" 

4600  ENTER  709  ;  G3 
4610  OUTPUT  709  i,,OPN" 

46Z0  G4-1 /(  S( 4, 1 )+S< 4 ,Z )*LOG( G3)+S(4 , 3)»L0G( G3>* 31-273. IS 
4630  PRINT  USING  Z10  t  G4 
4640  URIT  GZ  «  NEXT  I 
4650  GOTO  350 

4660  IF  GOS  THEN  GOTO  4770 
4670  FOR  1-1  TO  G1 
4680  OUTPUT  709  rCLSO” 

4690  OUTPUT  709  j"LS6" 

4700  OUTPUT  709  5 "F3N5RROR5T3" 

4710  ENTER  709  ;  G3 
47Z0  OUTPUT  709  5MOPN" 

4730  G4-1/(  S(  5, 1 >+S( 5 , Z ) *LOG< G3)+S(  5 , 3 >*LOG(  G3>“ 3)-Z73. 15 
4740  PRINT  USING  Z10  «  G4 
4750  URIT  GZ  Q  NEXT  I 
4760  GOTO  350 

4770  IF  GOG  THEN  GOTO  4880 
4780  FOR  1-1  TO  G1 
4790  OUTPUT  709  ;"CLS0" 

4800  OUTPUT  709  ^‘LS7,, 

4810  OUTPUT  709  1 "F3N5RROR5T3" 

48ZO  ENTER  709  i  G3 
4830  OUTPUT  709  r’OPN" 

4840  G4-1/(S(  16,1 >+S< 1 6 ,Z >*LOG( G3>+S( 16,3)*L0G(G3>* 3J-Z73. 15 
4850  PRINT  USING  Z10  1  G4 
4860  URIT  GZ  0  NEXT  I 
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4870  GOTO  350 

4880  IF  G<>7  THEN  GOTO  4990 
4890  FOR  1-1  TO  G1 
4900  OUTPUT  709  ;"CLS0” 

4910  OUTPUT  709  ;,,LS8“ 

4920  OUTPUT  709  ; "F3NSRft0R5T3“ 

4930  ENTER  709  ;  G3 
4940  OUTPUT  709  rOPN" 

4950  G4«1/(S< 17,1 >+S< 17, 2 )*LOG< G3>+S< 17,3)*L06<G3)A3)-Z73. 15 
4960  PRINT  USING  210  ?  G4 
4970  URIT  GZ  @  NEXT  I 
4980  GOTO  350 

4990  IF  G<>8  THEN  GOTO  5100 
5000  FOR  1-1  TO  G1 
5010  OUTPUT  709  5“CLS0" 

50Z0  OUTPUT  709  ;"LS9‘* 

5030  OUTPUT  709  ; "F3N5RA0R5T3" 

5040  ENTER  709  ;  G3 
5050  OUTPUT  709  rOPN" 

5060  G4-1 /( S( 6, 1 )+S(  6,2  >*LQG( G3)+S<  6, 3)*L0G(  G3) A 3)-Z73. 15 
5070  PRINT  USING  210  j  G4 
5080  UftIT  GZ  0  NEXT  I 
5090  GOTO  350 

5100  IF  G<>9  THEN  GOTO  5210 
5110  FOR  1-1  TO  G1 
5120  OUTPUT  709  i“CLSr 
5130  OUTPUT  709  I-LS2" 

5140  OUTPUT  709  ;  ,,F3N5RA0R5T3“ 

5150  ENTER  709  t  G3 
5160  OUTPUT  709  .•"OPN” 

5170  G4-1/( S( 7, 1 )+S( 7.2  >*LOG(  G3)+S(7, 3)*LOG( G3>* 3 )-Z73. 15 
5180  PRINT  USING  210  »  G4 
5190  WAIT  G2  0  NEXT  I 
5200  GOTO  350 

5210  IF  G0 10  THEN  GOTO  5320 
5220  FOR  1-1  TO  G1 
5Z30  OUTPUT  709  .-"CLSI" 

5240  OUTPUT  709  >"LS3M 

5Z50  OUTPUT  709  \ "F3NSRA0RST3" 

5260  ENTER  709  <  G3 
5270  OUTPUT  709  »"OPN“ 

5280  G4-1/( S( 8, 1 )+S(  8 , 2  >*LOG( G3)+S(  8 , 3>»LOG( G3) A 3 )-Z73. 15 
5290  PRINT  USING  210  <  64 
5300  UAIT  GZ  C  NEXT  I 
5310  GOTO  350 

5320  IF  G0 11  THEN  GOTO  5440 
5330  OUTPUT  709  rCLSr 
5340  FOR  1-1  TO  G1 
5350  OUTPUT  709  i"CLSr 
5360  OUTPUT  709  |"LS4" 
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5370  OUTPUT  709  ;  "F3N5RA0R5T3" 

5380  ENTER  709  j  63 
5390  OUTPUT  709  i"OPN" 

5400  G4-  1/(  S(  9, 1  >+S(  9 , Z  >*LOG(  63  >+S(  9 , 3)»LOG(  G3)*  3)-273.  15 
5410  PRINT  USING  Z10  j  G4 
54Z0  WRIT  GZ  8  NEXT  I 
5430  GOTO  350 

5440  IF  G0 12  THEN  GOTO  5550 
5450  FOR  1-1  TO  G1 
5460  OUTPUT  709  ;"CLSr 
5470  OUTPUT  709  j"LS5" 

5480  OUTPUT  709  j "F3N5RA0R5T3" 

5490  ENTER  709  !  G3 
5500  OUTPUT  709  s"OPN“ 

5510  G4-1/(S( 10, 1 >+S< 10, Z )*LOG( G3 >+S(  10, 3>*LOG(G3>* 3>-Z73. 15 
5520  PRINT  USING  Z10  i  G4 
5530  WRIT  GZ  8  NEXT  I 
5540  GOTO  350 

5550  IF  G0 13  THEN  GOTO  SEG0 
5560  FOR  1-1  TO  G1 
5570  OUTPUT  709  j"CLS1" 

5580  OUTPUT  709  i"LS6" 

5590  OUTPUT  709  } "F3N5RR0R5T3" 

5600  ENTER  709  i  G3 
5610  OUTPUT  709  «"OPN" 

5620  G4-1/(S( 1 1 ,1 )+S( 1 1 ,Z)*L0G(G3>+S( 1 1 ,3)*L0G(G3)*3)-Z73. 15 
5630  PRINT  USING  Z10  j  G4 
5640  WRIT  GZ  8  NEXT  I 
5650  GOTO  350 

5660  IF  G0 14  THEN  GOTO  5770 
5670  FOR  1-1  TO  G1 
5680  OUTPUT  709  {"CL SI" 

5690  OUTPUT  709  !MLS7" 

5700  OUTPUT  709  j "F3N5RA0R5T3" 

5710  ENTER  709  i  G3 
57Z0  OUTPUT  709  r'OPN" 

5730  G4*l/( S( 1 Z , 1 )+S< 1Z,Z)*L0G(G3>+S( 1Z , 3 >*LOG< G3 >A 3)-Z73. 15 
5740  PRINT  USING  Z10  j  G4 
5750  URIT  GZ  8  NEXT  I 
5760  GOTO  350 

5770  IF  G0 15  THEN  GOTO  5870 
5780  FOR  1-1  TO  G1 
5790  OUTPUT  709  j"CLS1" 

5800  OUTPUT  709  j"LS8" 

5810  OUTPUT  709  j "F3N5RA0R5T3" 

58Z0  ENTER  709  i  G3 
5830  OUTPUT  709  »"OPN" 

5840  G4-1/(S(  13,1  >+S<  1 3,Z  >*LOG<  G3HS(  1  3, 3  )*LOG<  G3  >*3>-Z73.  15 
5850  PRINT  USING  Z10  i  G4 
5860  URIT  GZ  8  NEXT  I 
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5870 
5880 
5890 
5900 
5910 
5920 
5930 
5940 
5950 
59  G0 
5970 
5980 
5990 
6000 
6010 
6020 
6030 
6040 
6050 
6060 
6070 
6080 
6090 
6100 
6110 
6120 
6130 
6140 
6150 
6160 
6170 
6180 
6190 
6200 
6210 
6220 
6230 
6235 
6240 
6250 
6260 
6270 
6280 
6290 
6300 
6310 
6320 
6330 
6340 
6350 


IF  G0 16  THEN  GOTO  5970 
FOR  1-1  TO  G1 
OUTPUT  709  i "CLS1" 

OUTPUT  709  ; “LS9" 

OUTPUT  709  i "F3N5RA0R5T3" 

ENTER  709  ;  G3 
OUTPUT  709  i "OPN" 

G4-1/(S< 14,1 )+S< 14,2 )*L0G( G3)+S< 14,3)*L0G(G3>*3)-273. 15 
PRINT  USING  210  ;  64 
UAIT  62  6  NEXT  I 
DISP  "NOSUCH  CHANNEL" 

GOTO  350 
I 

DATA  .0013217799,. 00025777 108,1. 7435609E- 7 
DATA  . 001 3546937 , . 0002530644 . 1 . 8275039E-7 
DATA  .0013105286, .00025873507, l .7603274E-7 
DATA  .001 348592 . . 00025330525 . 1 . 5758 1 63E-7 
DATA  .001 302763 . . 000261 17401,1. 6079Z75E-7 
DATA  . 001 31 81 985. . 000261 39129,1. 3640022E-7 
DATA  . 001 3343537,. 00025692 141 . 1 . G000615E-7 
DATA  . 00 1 3037492 . . 0002 5949 1 84 , 1 . 66634 1 E-7 
DATA  .001 3348872 , .00025676909 , 1 . 6407671 E-7 
DATA  .001 3526671 . . 00025602699 , I . 5699062E-7 
DATA  .0013006584,-00025973931 , 1 .6871901E-7 
DATA  .0013503701 ,.00025502833. 1.51 G8687E-7 
DATA  .0013124103,. 0002594037 . 1 . 701 2 86 1 E-7 
DATA  .0013027286, .00025987376. 1 .6883065E-7 
DATA  . 001 3343537. . 00025692 141,1. 600061 5E-7 
DATA  .0013545896, .00025492581 , 1 .5161931E-7 
DATA  .0013122181 . .00025971489 , 1 . 6281 56E-7 
DATA  .0013348872,-00025676909,1 .640767 IE-7 
DATA  . 001310042,. 00025963523, J.6717593E-7 
DATA  .001321 7799 , . 00025777 108,1. 7435609E-7 
I 

PRINTER  IS  707,80 
OUTPUT  707 
OUTPUT  707 

PRINT  "Teat  no.  ";TS 
Z-IP<  T ) 

Z1-IP(  ( T-Z)*60)  !  MIN 
IMAGE  "TIME", IX, ZD. ZD 
PRINT  USING  6270  ;  Z,Z1 
PRINT  Dt 

IMAGE  3D,"  gal/nin" 

PRINT  "FLOW  RATE:";®  PRINT  USING  6300  ;  FI/. 00006 30902 
IMAGE  3D. ZD."  Deg  C" 

PRINT  "BULK  UATER  TEMP";®  PRINT  USING  6320  ;  P4<3> 

PRINT  "PRANDTL  NO.  "; 

PRINT  USIN6  210  ;  B7 
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G360  PRINT 

6370  IMAGE  14X," STATION  1“ , 3X," STATION  Z"  , 3X, "STATION  3" , 3X." STATION  4" 
6380  PRINT  USING  6370 
6390  PRINT 

6400  PRINT  "PLATE  TEMP" 

6410  PRINT  “  6  0  cm  "5 

6420  IMAGE  5X. 3D.2D, GX, 3D. 2D, GX . 30. 2D, 6X, 3D. ZD 
6430  PRINT  USING  6420  ;  T< 1 ) ,T< 2 ) ,T( 3> ,T<4> 

6440  PRINT  "  6  .325  cxi";®  PRINT  USING  6420  !  R( 2 ) ,R( 4 > ,R(  6) ,R< 8 ) 

6450  PRINT  "  6  1.59  cm";®  PRINT  USING  6420  5  R(  1 > ,R(  3) ,R(  5) ,R( 7) 

6455  PRINT  "8  1.914  cm";®  PRINT  USING  6420  ;  TG( 1 > . T6< 2 > , T6( 3) , T6( 4) 
6460  IMAGE  8X,6D.ZD,5X,GD.ZD, 5X,6D. 2D.5X.6D.ZD 
6470  PRINT 

6480  PRINT  "HEAT  FLOW" 

6490  PRINT  "METER" 5®  PRINT  USING  646©  ;  H( 1 ) , H< 2 ) ,H( 3) ,H( 4) 

6500  PRINT  "X  .55";®  PRINT  USING  6460  ;  HI ( 1 > ,H1 ( 2 ) ,H1 ( 3) .HI ( 4 ) 

6520  PRINT 

6530  IMAGE  5X,70,6X,7D,6X,70,6X,7D 
6540  PRINT  "REYNOLDS  NUMBER" 

6550  PRINT  "  HYDR  DIAM";®  PRINT  USING  6530  ;  R5( 1 > ,RS(  2 ) ,R5(  3) ,R5( 4) 

6560  PRINT  "  ENTR  LEN  "i®  PRINT  USIN6  6530  i  R6< 1 > ,R6( 2 ) ,R6( 3 > ,R6( 4 ) 

6570  IMAGE  1 1X.6D.8X,GDf8X,60.8X,6D 
6580  PRINT 

6590  PRINT  “NUSSELT  NUMBER,  HYDR.  DIAM" 

6600  PRINT  "METER" 1 ®  PRINT  USING  6570  ;  N1 ( 1 ) ,N1 ( 2 ) ,N1 ( 3) ,N1 ( 4 ) 

6610  PRINT  "X  .55"}®  PRINT  USING  6570  5  NZ( 1 ) ,NZ< 2 > ,NZ( 3) ,NZ< 4 ) 

6630  PRINT 

6640  PRINT  "NUSSELT  NO.,  ENTR.  LEN." 

6650  PRINT  “  METER";®  PRINT  USING  6570  ;  N4< 1 ) ,N4(Z) ,N4( 3) ,N4(4) 

6660  PRINT  "  X  .55";®  PRINT  USING  6570  ;  N5( 1 > ,N5( 2 > ,N5( 3) ,N5< 4 > 

6680  RETURN 
6690  END 

6700  I  EXTERNAL  SUBROUTINE 
6710  !  TO  SET  UP  MOTOR  CONTROL 
6720  !  MOTOR  CONTROL  ADDRESS 
6730  I  IS  10 
6740  ! 

6750  DISP  "SET  PROBE  TO  (0,0)" 

6760  DISP  "PRESS  CONT  WHEN  READY"  e  PAUSE 

6770  !  SETUP  HP-IL/RSZ3ZC  INTERFACE 

678fl»  CONTROL  10,3  ;  6 

6790  CONTROL  10,4  ;  10 

6800  CONTROL  10,12  ;  94 

6810  CONTROL  10,13  ;  13 

6820  CONTROL  10,14  ;  63 

6830  CONTROL  10,16  ;  1 

6840  CONTROL  10,17  ;  13 

6850  CONTROL  10,9  ;  136 
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6880  I  UNITS  WILL  BE  CENTIMETERS 
6870  oi sp  "6650" 

6880  OUTPUT  10  5 “E" 

6890  DISP  "6670" 

6900  UAIT  2000 
6920  I 

6930  OUTPUT  10  * "0  12  — I2:IFA2<>P2  THEN  A2— A2" 

6940  WAIT  2000 

7190  OUTPUT  10  ; "200  P-1: INPUT  X,S6,S7" 

7200  OUTPUT  10  ; “205  60SUB  1000:A1-X:@:A2-0:@" 

7220  DISP  ''6970" 

7230  OUTPUT  10  ;  “2 10  GOSUB  500: STOP: FOR  Q-1  TO  <S7-1>" 

7240  WAIT  1000 

7250  !  OUTPUT  10  s"220  IF(P1-.5><-19.8  THEN240" 

7260  01 SP  "70 10“ 

7280  OUTPUT  10  s"230  I2-S6: 6:G0SUB  500:STOP:NEXT  Q" 

7300  OUTPUT  10  i "240  A1 -0: 6: A2-0: @:END" 

7310  DISP  "7060" 

7595  OUTPUT  10  j"500  P6-INT(P1  •  10*2+ . 5 >/ 10*2 :P7-INT< P2*  10*2+ . 5 )/ lO^" 

7596  UAIT  2000 

7600  OUTPUT  10  i "510  PRINT  P6:PRINT  P7: RETURN" 

7620  DISP  "7370“ 

7630  OUTPUT  10  t " 1000  C 1 - .000635: C2-. 000635 : R1 - 1 : R2-R1 " 

7650  DISP  "7400" 

7660  OUTPUT  10  1010  V 1 -900: V2-V 1  :D1- 10:  RETURN" 

7670  01 SP  "7420" 

7680  RETURN 
7690  END 
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APPENDIX  B:  TYPICAL  DATA  OUTPUT,  TEST  066 


Test  no.  0065 
TIME  13:47 
6  BUG  86 

FLOW  RATE: 400  gal/nin 
BULK  WATER  TEMP  8.Z8  Deg  C 
PRANOTL  NO.  10. M 


STATION  1 

STATION  Z 

STATION  3  STATION  4 

PLATE  TEMP 

1  0  Cfl 

3.73 

4.93 

4.00 

5.  1 1 

6  . 3Z5  cn 

6.19 

6.6Z 

6.47 

6.59 

6  1.59  cn 

6.53 

6.87 

6.77 

6.86 

e  1.914  cn 

6.6Z 

6.93 

6.84 

6.93 

HEAT  FLOW 

METER 

-364.97 

-198.50 

-Z48.68 

-181.88 

X  .55 

-3634.36 

-2657.30 

-3Z10.97 

-Z9Z8.S9 

REYNOLDS  NUMBER 

HYOR  01 AM  1 164 II 

1 16890 

116760 

1 16830 

ENTR  LEN 

3995518 

41Z9839 

4Z40749 

4357679 

NUSSELT  NUMBER,  HYDR.  DIAM 
METER  1 1 1 

76 

89 

7Z 

X  .55 

1109 

1016 

1150 

1155 

NUSSELT  NO. , 
METER 

ENTR.  LEN. 
38Z1 

Z68Z 

3Z36 

Z676 

X  .55 

38053 

35898 

41784 

43094 

TIME  13.7897958333 
AIR  TEMP  Z5.73 
FLUME  ENTRANCE  TEMP  8.68 
FLUME  EXIT  TEMP  8.59 
COLD  GLYCOL  ENTR  JZ7. 1 1 
COLD  GLYCOL  EXIT  -Z5.47 
HOT  GLYCOL  ENTR  11.54 
HOT  GLYCOL  EXIT  11. Z5 

FLUME  ENERGY  -9.Z65 
COLD  GYLCOL  EN  I1.47Z 
HOT  GYLCOL  EN  0.000 
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STATION 


XCOORD 

YCOORD 

0.00 

.  10 

0.00 

1  .  10 

0.00 

2.10 

0.00 

3.  10 

0.00 

4.  10 

0.00 

5.  10 

0.00 

6.  10 

0.00 

7.  10 

0.00 

8.  10 

0.00 

9.10 

0.00 

10.  t0 

0.00 

11.10 

0.00 

12.10 

0.00 

13.  10 

TEMP(C) 

VEL( M/S 

8.  19 

.02 

8.33 

.02 

8.35 

.02 

8.42 

.02 

8.40 

.02 

8.42 

.02 

8.43 

.02 

8.45 

.02 

8.46 

.02 

8.47 

.02 

8.48 

.02 

8.48 

.02 

8.49 

.02 

8.50 

.02 

K  OF  WATER  .579704283168 
HEAT  TRAN  CO  2090.96338324 
HYDR  01 AM  .307378009109 
VISCOSITY  OF  WATER  1 .4104400871 8E-6 
VELOCITY  OF  WATER  .534164761663 
DENSITY  OF  WATER  999.556769831 


STATION 

XCOORD 

YCOORD 

TEMPI C ) 

VEL( M/S 

2 

0.00 

.  10 

8. 17 

.02 

2 

0.00 

1 .  10 

8.25 

.02 

2 

0.00 

2.10 

8.31 

.02 

2 

0.00 

3.  10 

8.37 

.02 

2 

0.00 

4.  T0 

8.39 

.02 

2 

0.00 

5.  10 

8.41 

.02 

2 

0.00 

6.  10 

8.40 

.02 

2 

0.00 

7. 10 

8.43 

.02 

2 

0.00 

8.  10 

8.44 

.02 

2 

0.00 

9.10 

8.47 

.02 

2 

0.00 

10.  10 

8.47 

.02 

2 

0.00 

11.10 

8.47 

.02 

2 

0.00 

12.10 

8.48 

.02 

2 

0.00 

13.  10 

8.49 

.02 

K  OF  WATER  .579951897005 
HEAT  TRAN  CO  1917.05079022 
HYDR  DIAM  .307378009109 
VISCOSITY  OF  WATER  I . 40466247544E-6 
VELOCITY  OF  WATER  .534164761663 
DENSITY  OF  WATER  999.54463703 
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STATION 

XCOORD 

Y COORD 

TEMP( C ) 

VELU1/S 

3 

0.00 

.  10 

8.25 

.02 

3 

0.00 

1 .  10 

8.25 

.02 

3 

0.ee 

2.10 

8.32 

.02 

3 

0.00 

3.  10 

8.34 

.02 

3 

0.00 

4.  10 

8.33 

.02 

3 

0.00 

5.  10 

8.37 

.02 

3 

0.00 

6.  10 

8.39 

.02 

3 

0.00 

7.  10 

8.41 

.02 

3 

0.00 

8.  10 

8.40 

.02 

3 

0.00 

9.10 

8.42 

.02 

3 

0.00 

10. 10 

8.44 

.02 

3 

0.00 

11.10 

8.45 

.02 

3 

0.00 

12.  10 

8.47 

.02 

3 

0.00 

13.  10 

8.44 

.02 

K  OF 

WATER  .579885136096 

HEAT 

TRAN  CO  2170. 

35066979 

HYDR 

DIAM  .307378009109 

VISCOSITY  OF  WATER 

!  1 .4062 173382 3E- 6 

VELOCITY  OF  WATER 

.534164761663 

DENSITY  OF  WATER 

999.547930213 

STATION  XCOORD 

YCOORD 

TEMPI C ) 

VEL( M/S 

4 

0.00 

.  10 

8.  19 

.02 

4 

0.00 

1.10 

8.23 

.02 

4 

0.00 

2.  10 

8.28 

.02 

4 

0.00 

3. 10 

8.32 

.02 

4 

0.00 

4. 10 

8.35 

.02 

4 

0.00 

5.  10 

8.32 

.02 

4 

0.00 

6.  10 

8.37 

.02 

4 

0.00 

7.  10 

8.35 

.02 

4 

0.00 

8. 10 

8.38 

.02 

4 

0.00 

9.10 

8.43 

.02 

4 

0.00 

10. 10 

8.40 

.02 

4 

0.00 

11.10 

8.43 

.02 

4 

0.00 

12.  10 

8.43 

.02 

4 

0.00 

13.  10 

8.45 

.02 

K  OF  WATER 

!  .579921040138 

HEAT  TRAN 

CO  2179.75455364 

HYDR  DIAM 

.307378009109 

VISCOSITY  OF  WATER  1 . 40538086839E-6 
VELOCITY  OF  WATER  .534164781683 
DENSITY  OF  WATER  999.546161154 
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APPENDIX  C:  TEST  DATA  AND  CALCULATED  PARAMETERS 


Notes: 

(1)  Lower  thermistor  in  flume  bed. 

(2)  Upper  thermistor  in  flume  bed. 

(3)  Calculated  bed  surface  temperature. 

(4)  Reynolds  number  using  velocity  correction. 

(5)  Velocity  correction  is  not  applicable  to  low  flow  rates. 

(6)  Entries  116-119  were  not  included;  these  data  were  not  at  steady  state. 


Fits*  bad  ta^aratura  abova  4  dagraaa  C,  turbulant  ftoa 


Entry  Station 
No.  No. 

Taaparatura  at 
(1)  (2)  (3) 

(daaraaa  C) 

Flan 

tats 

(GPU) 

ttotar  bulk 
Oapth  Taap 
(ca)  <C> 

Pr 

no 

Hydraulic  Dlaau 
ta  Na<4) 

t'tf 

Mu 

Entranca  Length 

Ha  Nu 

Haat 

Flux 

(U/a*2) 

1 

2 

6.50 

6.83 

6.92 

170 

14.5 

10.7 

9.1 

56362 

75999 

489 

1.956*06 

1.808*04 

3608.6 

2 

3 

6.27 

6.54 

6.61 

170 

14.5 

10.7 

9.1 

56994 

75975 

364 

2.000*06 

1.308*04 

2933.9 

3 

4 

6.34 

6.67 

6.75 

170 

14.5 

10.6 

9.1 

56329 

75909 

470 

2.058*06 

1.821-04 

3594.3 

4 

2 

6.58 

6.91 

7.00 

170 

14.5 

10.7 

9.1 

56394 

76054 

497 

1.958*06 

1.836*04 

3608.6 

3 

3 

6.35 

6.64 

6.71 

170 

14.5 

10.7 

9.1 

56342 

75948 

405 

2.008*06 

1.53E-04 

3151.2 

« 

4 

6.42 

6.76 

6.85 

170 

14.5 

10.7 

9.1 

56378 

75999 

492 

2.068*06 

1.91E-04 

3703.2 

7 

2 

6.33 

6.65 

6.73 

170 

14.0 

10.5 

9.2 

57018 

77051 

466 

2.018*06 

1.758*04 

3499.2 

8 

3 

6.05 

6.37 

6.45 

170 

14.0 

10.5 

9.2 

57010 

76981 

431 

2.058*06 

1.668*04 

3477.2 

9 

4 

6.06 

6.40 

6.49 

170 

14.0 

10.4 

9.2 

56961 

76915 

467 

2.11E«06 

1.85E*04 

3703.2 

10 

2 

6.28 

6.60 

6.68 

170 

14.0 

10.4 

9.2 

56927 

76906 

447 

2.008*06 

1.75E*04 

3499.2 

11 

3 

6.00 

6.32 

6.40 

170 

14.0 

10.4 

9.2 

56894 

76800 

434 

2.058*06 

1.678*04 

3477.2 

12 

4 

6.03 

6.38 

6.47 

170 

14.0 

10.4 

9.2 

56861 

76768 

486 

2.118*06 

1.92E*04 

3812.2 

13 

2 

7.08 

7.41 

7.50 

220 

13.5 

10.5 

9.2 

75149 

106649 

586 

2.738*06 

2.238*04 

3608.6 

14 

3 

6.82 

7.13 

7.23 

220 

13.5 

10.5 

9.2 

75149 

106580 

535 

2.798*06 

2.10E*04 

3585.8 

15 

4 

6.91 

7.27 

7.36 

220 

13.5 

10.5 

9.2 

75061 

106479 

617 

2.871*06 

2.488*04 

3921.1 

16 

2 

7.08 

7.41 

7.50 

220 

13.5 

10.5 

9.2 

75171 

106682 

584 

2.738*06 

2.228*04 

3608.6 

17 

3 

6.84 

7.16 

7.24 

220 

13.5 

10.5 

9.2 

75127 

1065a 

522 

2.798*06 

2.058*04 

3477.2 

18 

4 

6.91 

7.27 

7.36 

220 

13.5 

10.5 

9.2 

75105 

106544 

613 

2.878*06 

2.478*04 

3921.1 

19 

2 

4.28 

4.55 

4.62 

100 

14.5 

10.5 

9.2 

32952 

39922 

257 

1.108*06 

9.498*03 

2952.5 

20 

3 

4.03 

4.27 

4.33 

100 

14.5 

10.5 

9.2 

32901 

39925 

215 

1.138*06 

8.178*03 

2607.9 

21 

4 

3.93 

4.22 

4.29 

100 

14.5 

10.5 

9.2 

33010 

MAUI 

JrWV 

258 

1.168*06 

1.018*04 

3158.6 

22 

2 

4.29 

4.54 

4.60 

100 

16.5 

10.4 

9.2 

32943 

39907 

238 

1.108*06 

8.788*03 

2733.8 

23 

3 

4.02 

4.27 

4.35 

100 

14.5 

10.4 

9.2 

32943 

99874 

226 

1.136*06 

8.578*03 

2716.5 

24 

4 

3.91 

4.19 

4.26 

100 

14.5 

10.4 

9.2 

32876 

39776 

254 

1.168*06 

9.888*03 

3049.7 

25 

2 

4.36 

4.59 

4.65 

100 

14.5 

10.5 

9.2 

32952 

39925 

220 

1.108*06 

8.128*03 

2515.1 

26 

3 

4.06 

4.31 

4.37 

100 

14.5 

10.4 

9.2 

32933 

99866 

228 

1.138*06 

8.448*03 

2716.5 

27 

4 

3.95 

4.22 

4.29 

100 

14.5 

10.4 

9.2 

32909 

99818 

244 

1.168*06 

9.528*03 

2940.8 

28 

2 

6.58 

6.86 

6.93 

150 

14.0 

10.5 

9.2 

50518 

430 

1.708*06 

1.618*04 

3061.8 

29 

3 

6.35 

6.61 

6.68 

ISO 

14.0 

10.5 

9.2 

90SQ3 

44440 

371 

1.828*06 

1.438*04 

2825.2 

30 

4 

6.37 

6.66 

6.73 

150 

14.0 

10.4 

9.2 

50245 

66995 

426 

1.878*06 

1.688*04 

3158.6 

31 

2 

6.55 

6.84 

6.91 

ISO 

14.0 

10.4 

9.2 

50186 

46343 

as 

1.778*06 

1.708*04 

3171.2 

32 

3 

6.32 

6.59 

6.66 

ISO 

14.0 

10.3 

9.2 

50113 

66191 

990 

1.828*06 

1.538*04 

2933.9 

33 

4 

6.33 

6.64 

6.72 

ISO 

14.0 

10.3 

9.2 

50113 

44202 

468 

1.878*06 

1.848-04 

3376.5 

34 

2 

7.11 

7.41 

7.49 

200 

13.5 

10.4 

9.2 

68199 

95153 

sa 

2.488*06 

2.098*04 

3280.5 

35 

3 

6.89 

7.19 

7.27 

200 

13.5 

10.4 

9.2 

68099 

94949 

514 

2.538*06 

2.018*04 

3259.8 

36 

4 

6.93 

7.24 

7.32 

200 

13.5 

10.4 

9.2 

68099 

95052 

536 

2.611*06 

2.168*04 

3376.5 

37 

2 

7.12 

7.41 

7.48 

200 

13.5 

10.4 

9.2 

41099 

95092 

532 

2.478*06 

2.038*04 

3171.2 

38 

3 

6.90 

7.19 

7.26 

200 

13.5 

10.4 

9.2 

WmKrw 

94979 

05 

2.538*06 

1.948*04 

3151.2 

39 

4 

6.95 

7.26 

7.34 

200 

13.5 

10.4 

9.2 

40020 

94937 

5*7 

2.608*06 

2.208*04 

3376.5 

40 

2 

7.7* 

8.0i 

8.09 

247 

14.0 

10.3 

9.2 

12520 

119189 

680 

2.978*06 

2.546*04 

3061.8 

41 

3 

7.48 

7.85 

7.9* 

247 

13.9 

10.3 

9.2 

82779 

119573 

837 

3.078*06 

3.238*04 

4020.5 

42 

4 

7.48 

7.81 

7.89 

247 

13.8 

10.3 

9.2 

81016 

119952 

737 

3.178*06 

2.938*04 

3594.3 

43 

2 

7.78 

8.06 

8.13 

247 

14.0 

10.3 

9.2 

02520 

119201 

692 

2.978*06 

2.588*04 

3061.8 

44 

3 

7.52 

7.90 

8.00 

247 

13.9 

10.4 

9.2 

02027 

119664 

872 

3.078*06 

3.368*04 

4129.1 

45 

4 

7.52 

7.86 

7.95 

247 

13.8 

10.3 

9.2 

85040 

120005 

773 

3.178*06 

5.078*04 

3703.2 

46 

2 

8.89 

9.10 

9.15 

300 

13.8 

10.7 

9.1 

101918 

152522 

7a 

3.748*06 

2.008*04 

2296.4 

47 

3 

8.75 

9.00 

9.06 

300 

13.7 

10.7 

9.1 

102300 

153124 

823 

3.868*06 

3.191*04 

2716.5 

41 


tntry  Station  Twparaturt  ot 

No.  No.  (1)  _  (2)  (3) 

(cUsrooo  C) 

H« 

lata 

(EM) 

Watar  bulk 
Oaptk  Tap 
(ea)  <C) 

Pr 

no 

Hydraulic  Dlaaatar 
•a  la<4)  Nu 

Cntranca  lanath 

la  Nu 

Haat 

Flux 

C«/«-2J 

48 

4 

8.80 

9.0 

9.09 

300 

13.7 

9.1 

102270 

15300 

774 

3.971*06 

3.098*04 

2303.1 

49 

2 

8.94 

9.17 

9.0 

300 

13.8 

10.4 

9.1 

101712 

152213 

820 

3.748*04 

3.00*04 

2296.4 

SO 

3 

8.0 

9.07 

9.13 

300 

13.7 

9.1 

102152 

152908 

01 

3.878*04 

3.308*04 

2607.9 

SI 

4 

8.87 

9.10 

9.14 

300 

13.7 

10.4 

9.1 

102152 

152918 

80 

3.978*04 

3.38*04 

88.1 

S2 

2 

8.89 

9.18 

9.0 

350 

13.5 

10.5 

9.2 

119520 

183357 

1244 

4.458*04 

4.748*04 

3171.1 

S3 

3 

8.84 

9.11 

9.17 

350 

13.3 

10.5 

9.2 

119520 

18330 

1001 

4.578*04 

3.918*04 

2716.5 

54 

4 

8.92 

9.15 

9.21 

350 

13.5 

10.5 

9.2 

119489 

183291 

954 

4.708*04 

3.848*04 

88.1 

ss 

2 

9.01 

9.22 

9.27 

350 

13.5 

10.5 

9.2 

119555 

183622 

909 

4.458*04 

3.48*04 

2296.4 

so 

3 

8.89 

9.14 

9.20 

350 

13.5 

10.5 

9.2 

119520 

183334 

108 

4.578*04 

4.00*04 

2714.5 

sr 

4 

8.0 

9.17 

9.0 

350 

13.5 

10.5 

9.2 

119520 

183354 

1345 

4.708*04 

3.40*04 

3485.4 

so 

2 

9.14 

9.37 

9.42 

400 

13.3 

10.6 

9.2 

137751 

215860 

990 

5.188*04 

3.798*04 

2296.4 

59 

3 

9.05 

9.31 

9.38 

400 

13.3 

10.4 

9.2 

137751 

215837 

1147 

5.38*04 

4.598*04 

288.2 

00 

4 

9.10 

9.33 

9.39 

400 

13.3 

10.6 

9.1 

1377*1 

215911 

1030 

5.471*04 

4.208*04 

88.1 

01 

2 

9.22 

9.44 

9.50 

400 

13.3 

10.5 

9.2 

137711 

215827 

1119 

5.198*04 

4.298*04 

248.7 

02 

3 

9.12 

9.37 

9.43 

400 

13.3 

10.4 

9.2 

137751 

215844 

1100 

5.338*04 

4.648*04 

2716.5 

03 

4 

9.17 

9.40 

9. a 

400 

13.2 

10.4 

9.1 

138274 

214744 

1099 

3.518*04 

4.468*04 

238.1 

44 

2 

9.28 

9.49 

9.54 

450 

13.0 

10.4 

9.2 

155045 

248691 

1318 

5.948*04 

3.111*04 

2296.4 

45 

3 

9.19 

9.45 

9.52 

450 

13.0 

10.4 

9.2 

155754 

248481 

1404 

6.12M4 

6.398*04 

288.2 

64 

4 

9.24 

9.47 

9.53 

450 

12.9 

10.4 

9.2 

154353 

249913 

1422 

4.348*04 

5.848*04 

258.1 

42 

2 

9.27 

9.48 

9.53 

450 

13.0 

10.3 

9.2 

155443 

248334 

1347 

5.958*04 

5.308*04 

2296.3 

48 

3 

9.18 

9.43 

9.49 

450 

13.0 

10.4 

9.2 

155709 

248391 

1520 

4.128*04 

6.08*04 

2716.3 

49 

4 

9.23 

9.45 

9.51 

450 

12.9 

10.4 

9.2 

154490 

249735 

1200 

4.348*04 

5.268*04 

896.2 

232 

1 

5.50 

3.77 

5.0 

200 

31.4 

13.4 

8.4 

45681 

58308 

268 

1.098*04 

6.978*03 

2904.8 

233 

2 

4.42 

4.84 

6.92 

200 

31.4 

13.4 

8.4 

45451 

58941 

283 

1.108*04 

7.568*03 

284.4 

234 

3 

4.38 

4.58 

6.43 

200 

31.4 

13.4 

8.4 

45491 

58493 

224 

1.138*04 

6.168*03 

2173.2 

235 

4 

4.34 

6.59 

6.45 

200 

31.4 

13.4 

8.4 

45691 

282 

1.168*04 

7.958*03 

2723.0 

234 

1 

5.50 

5.0 

5.91 

400 

30.5 

9.4 

9.4 

83594 

125273 

670 

2.058*04 

1.78*04 

3550.3 

237 

2 

4.12 

6.34 

6.42 

400 

30.5 

9.4 

9.4 

83594 

125489 

575 

2.18*04 

1.548*04 

284.4 

238 

3 

5.90 

6.16 

6.22 

400 

30.5 

9.4 

9.4 

83594 

125404 

50 

2.188*04 

1.618*04 

288.2 

239 

4 

5.98 

6.24 

6.31 

400 

30.5 

9.4 

9.4 

83594 

125440 

599 

2.248*04 

1.708*04 

201.9 

240 

1 

4.35 

6.70 

6.78 

300 

15.0 

9.2 

9.4 

93434 

138204 

819 

3.18*04 

2.878*04 

378.5 

241 

2 

4.0 

7.07 

7.15 

300 

15.0 

9.2 

9.4 

93434 

13830 

448 

3.28*04 

2.418*04 

284.4 

242 

3 

4.68 

4.94 

7.01 

300 

15.0 

9.2 

9.6 

93434 

13820 

758 

3.38*04 

2.818*04 

3151.2 

243 

4 

4.78 

7.05 

7.12 

300 

15.0 

9.2 

9.4 

93434 

138321 

74* 

3.448*04 

2.88*04 

2940.8 

244 

1 

4.19 

4.53 

4.61 

400 

15.5 

8.5 

9.8 

120315 

185358 

1034 

4.018*04 

3.58*04 

387.9 

24S 

2 

4.42 

4.87 

6.93 

400 

15.5 

S.S 

9.8 

120915 

185502 

935 

4.18*04 

3.311*04 

2733.8 

244 

3 

4.47 

6.77 

6.0 

400 

15.5 

0.5 

9.8 

12015 

10442 

1094 

4.268*04 

3.88*04 

389.8 

247 

4 

4.59 

6.84 

4.93 

400 

15.5 

0.5 

9.8 

12015 

185500 

1002 

4.38*04 

3.78*04 

2940.8 

248 

1 

4.21 

4.41 

4.71 

400 

15.2 

0.5 

9.8 

121574 

187409 

185 

4.18*04 

4.38*04 

4308.4 

249 

2 

4.59 

4.84 

4.90 

400 

15.2 

0.5 

9.8 

121530 

187517 

90S 

4.28*04 

3.28*04 

2733.8 

230 

3 

4.42 

4.74 

4.0 

400 

15.1 

8.4 

9.8 

121792 

187902 

1113 

4.378*04 

4.118*04 

3477.2 

2S1 

4 

6.53 

4.82 

4.09 

400 

15.0 

0.5 

9.8 

122248 

10727 

1042 

4.538*04 

3.98*04 

3158.6 

252 

1 

5.41 

5.79 

5.0 

500 

14.9 

7.5 

10.1 

148824 

234499 

1329 

5.088*04 

4.68*04 

4088.3 

2S3 

2 

5.79 

4.0 

4.09 

500 

14.9 

7.5 

10.1 

148024 

234915 

983 

5.258*04 

3.58*04 

284.4 

2S4 

3 

5.44 

5.94 

4.04 

500 

14.8 

7.5 

10.1 

149321 

237420 

181 

5.438*04 

4.68*04 

3477.2 

2SS 

4 

5.74 

4.0 

4.10 

500 

14.7 

7.5 

10.1 

149020 

239491 

1099 

5.48*04 

4.28*04 

2940.8 

2S4 

1 

5.48 

5.0 

5.90 

500 

14.9 

7.4 

10.2 

148445 

234400 

1234 

5.00*04 

4.38*04 

387.9 

257 

2 

5.48 

5.0 

5.91 

500 

14.9 

7.4 

10.2 

148445 

tmiim 

1242 

5.238*04 

4.58*04 

3717.9 

238 

3 

5.49 

5.98 

4.05 

500 

14.8 

7.4 

10.2 

149140 

237314 

1177 

5.48*04 

4.38*04 

3151.2 

259 

4 

5.81 

4.07 

4.14 

500 

14.7 

7.4 

10.2 

149438 

230200 

1121 

3.48*04 

4.38*04 

201.9 

Huh  bad  ta^aratura  ateva  4  < 

laaraaa  C 

lan  flea  rataa 

214 

1 

5.05 

5.10 

5.11 

o 

29.0 

32.3 

5.1 

9321 

(5) 

13 

1.88*0 

3.48*02 

537.9 

217 

2 

3.99 

4.08 

4.10 

o 

29.0 

32.3 

5.1 

9321 

24 

1.908*0 

6.88*02 

984.2 

218 

3 

4.20 

4.0 

6.24 

o 

29.0 

32.3 

5.1 

9321 

9 

1.948*0 

2.338*02 

326.0 

219 

4 

4.44 

4.51 

4.0 

a 

29.0 

32.3 

5.1 

9321 

19 

2.08*0 

5.658*02 

78.4 

220 

1 

5.04 

5.08 

3.09 

45 

31.0 

32.0 

5.1 

15942 

11 

3.048*0 

2.88*02 

430.3 

221 

2 

4.04 

4.13 

4.15 

43 

31.0 

32.0 

5.1 

15942 

8 

3.198*0 

4.98*02 

984.2 

222 

3 

4.24 

4.28 

4.29 

45 

31.0 

32.0 

5.1 

15942 

4 

3.298*0 

1.58*02 

217.3 

223 

4 

4.43 

4.50 

4.52 

45 

31.0 

32.0 

5.1 

15942 

20 

3.398*0 

5.78*02 

78.4 

224 

1 

7.53 

7.45 

7.68 

a 

30.3 

31.7 

5.1 

2404 

33 

5.378*0 

9.448*02 

1291.0 

223 

2 

8.0 

0.0 

8.34 

75 

30.3 

31.7 

5.1 

2404 

31 

3.50*0 

8.48*02 

1093.5 

42 


Entry  Station 
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9r 
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Entrance  length 

Meat 

No. 

NO. 

(2) 

(3) 

(at# 

Oaptft  Tmp 

na 

la 

la<4> 

Hu 

■a 

Du 

Mum 

(daireee  C) 

(099) 

(ON) 

<c> 

(tl/a*2) 

226 

3 

8.26 

8.29 

8.30 

73 

30.3 

31.7 

3.1 

2606 

9 

3.731*03 

2.590*02 

326.0 

227 

6 

8.36 

8.61 

8.63 

73 

30.3 

31.7 

3.1 

2606 

21 

3.891*05 

6.20*02 

762.4 

228 

1 

3.28 

3.48 

3.53 

100 

30.6 

32.9 

3.0 

36443 

51 

7.011*03 

1.30*0 

2151.7 

229 

2 

6.17 

6.31 

6.33 

100 

30.6 

32.9 

3.0 

36443 

38 

7.290*03 

1.00*0 

1330.9 

230 

3 

6.01 

6.07 

6.09 

100 

30.6 

32.9 

3.0 

36443 

16 

7.471*05 

4.510*02 

652.0 

231 

6 

3.86 

3.98 

6.01 

100 

30.6 

32.9 

5.0 

36443 

32 

7.670*03 

9.270*02 

1307.0 

Muaa  bad  t 

aneratura  balou  4  detract  c 

turbulent  flow 

70 

2 

1.78 

2.03 

2.12 

200 

13.3 

4.8 

11.1 

57333 

78038 

341 

2.090*06 

2.00*06 

2932.3 

71 

3 

1.33 

1.78 

1.86 

200 

13.3 

4.8 

11.1 

37333 

77997 

415 

2.140*06 

1.620*06 

2499.2 

72 

6 

1.37 

1.85 

1.92 

200 

13.3 

4.8 

11.1 

37497 

77961 

324 

2.200*06 

2.110*06 

309.7 

73 

2 

1.73 

2.00 

2.07 

200 

13.3 

4.8 

11.1 

37479 

77967 

337 

2.090*06 

2.00*0 

2932.3 

7* 

3 

1.30 

1.74 

1.80 

200 

13.3 

4.8 

11.1 

57533 

77989 

627 

2.140*06 

1.670*0 

2607.9 

73 

6 

1.36 

1.81 

1.88 

200 

13.3 

4.9 

11.1 

37531 

78032 

493 

2.200*06 

1.90*0 

200.8 

76 

1 

1.12 

1.47 

1.33 

100 

15.0 

7.8 

10.0 

29947 

33636 

SIS 

9.490*05 

1.110*0 

3765.3 

77 

2 

2.06 

2.30 

2.36 

100 

13.0 

7.8 

10.0 

29947 

33349 

252 

9.890*03 

9.130*0 

2624.4 

78 

3 

1.76 

1.98 

2.06 

100 

13.0 

7.8 

10.0 

aaau  m 

CrW 

33536 

236 

1.010*06 

8.70*0 

2607.9 

79 

6 

1.53 

1.81 

1.88 

100 

13.0 

7.8 

10.1 

29920 

33438 

231 

1.040*06 

9.600*0 

2831.9 

80 

1 

1.02 

1.40 

1.49 

100 

13.0 

7.8 

10.0 

29938 

33437 

339 

9.40*0 

1.10*0 

600.3 

81 

2 

1.96 

2.20 

2.26 

100 

13.0 

7.8 

10.0 

29974 

35373 

266 

9.80*0 

8.90*0 

2624.4 

82 

3 

1.70 

1.96 

2.03 

100 

15.0 

7.8 

10.1 

29911 

33663 

237 

1.010*0 

9.50*0 

2825.2 

as 

6 

1.30 

1.77 

1.86 

100 

13.0 

7.8 

10.0 

29929 

33463 

238 

1.00*0 

9.80*0 

2960.8 

86 

1 

1.72 

1.97 

2.03 

130 

14.3 

6.3 

10.6 

43627 

36011 

323 

1.40*0 

1.10*0 

2689.6 

83 

2 

2.63 

2.70 

2.76 

130 

14.5 

6.3 

10.6 

43614 

36113 

600 

1.310*0 

1.40*0 

2733.8 

86 

3 

2.28 

2.47 

2.32 

130 

14.3 

6.3 

10.6 

36033 

20 

1.30*0 

1.00*0 

200.6 

87 

6 

2.26 

2.30 

2.36 

130 

14.3 

6.3 

10.6 

43376 

36024 

30 

1.390*0 

1.60*0 

2614.0 

88 

1 

1.66 

1.90 

1.96 

ISO 

14.3 

6.3 

10.6 

43574 

35924 

336 

1.40*0 

1.20*0 

2797.2 

8V 

2 

2.33 

2.61 

2.68 

ISO 

14.3 

6.2 

10.6 

43333 

53986 

413 

1.308*0 

1.30*0 

2863.1 

90 

2.19 

2.60 

2.43 

130 

14.3 

6.2 

10.6 

43347 

aaiua 

311 

1.340*0 

1.10*0 

201.9 

91 

6 

2.18 

2.66 

2.31 

130 

14.3 

6.2 

10.6 

43360 

53996 

390 

1.390*0 

1.510*0 

201.9 

92 

1 

1.96 

2.33 

2.62 

230 

14.0 

3.1 

11.0 

71199 

100196 

004 

2.670*0 

2.90*0 

4193.8 

93 

2 

2.36 

2.79 

2.83 

230 

14.0 

3.1 

11.0 

71177 

100278 

629 

2.360*0 

2.30*0 

2733.8 

96 

3 

2.36 

2.66 

2.74 

230 

14.0 

3.1 

11.0 

71133 

100813 

764 

2.60*0 

2.930*0 

3477.2 

95 

6 

2.38 

2.67 

2.76 

230 

14.0 

3.0 

11.1 

71088 

100111 

704 

2.60*0 

2.70*0 

3158.6 

96 

1 

1.91 

2.26 

2.36 

230 

14.0 

3.0 

11.1 

71063 

99978 

717 

2.660*0 

2.60*0 

370.3 

97 

2 

2.30 

2.76 

2.80 

230 

14.0 

3.0 

11.1 

71043 

100039 

606 

2.350*0 

2.260*0 

2626.4 

98 

3 

2.29 

2.60 

2.68 

230 

14.0 

5.0 

11.1 

70999 

99937 

743 

2.620*0 

2.80*0 

330.5 

99 

6 

2.36 

2.62 

2.69 

230 

14.0 

3.0 

11.1 

70934 

99893 

60 

2.690*0 

2.60*0 

309.7 

100 

1 

2.41 

2.76 

2.86 

300 

13.7 

3.0 

11.1 

86264 

123399 

80 

3.00*0 

3.10*0 

370.5 

101 

2 

2.90 

3.13 

3.21 

300 

13.7 

3.1 

11.0 

86298 

123604 

70 

3.130*0 

2.810*0 

2733.8 

102 

3 

2.70 

3.01 

3.09 

300 

13.7 

3.1 

11.0 

86379 

123691 

866 

3.20*0 

3.20*0 

330.5 

103 

6 

2.73 

3.03 

3.10 

300 

13.6 

5.1 

11.0 

86760 

12018 

753 

3.350*0 

3.00*0 

309.7 

106 

1 

2.38 

2.96 

3.03 

300 

13.7 

3.2 

11.0 

86869 

123970 

910 

3.060*0 

3.30*0 

3873. 1 

103 

2 

3.07 

3.33 

3.60 

300 

13.7 

3.2 

11.0 

866ZS 

126173 

on 

3.170*0 

3.00*0 

280.1 

106 

3 

2.86 

3.18 

3.26 

300 

13.7 

3.2 

11.0 

86677 

12016 

099 

3.20*0 

3.40*0 

3477.2 

107 

6 

2.91 

3.20 

3.27 

300 

13.6 

3.2 

11.0 

87060 

126866 

806 

3.370*0 

3.20*0 

3138.6 

108 

1 

0.63 

1.06 

1.16 

120 

13.0 

6.4 

10.5 

36410 

42028 

40 

1.110*0 

1.30*0 

4411.0 

109 

2 

1.68 

1.71 

1.77 

120 

13.0 

6.3 

10.6 

34378 

62066 

20 

1.10*0 

1.00*0 

2315.1 

110 

3 

1.22 

1.31 

1.30 

120 

13.0 

6.3 

10.6 

34314 

41936 

01 

1.100*0 

1.310*0 

3131.2 

111 

6 

1.17 

1.63 

1.32 

120 

13.0 

6.3 

10.6 

4190 

336 

1.210*0 

1.20*0 

309.7 

112 

1 

0.33 

0.92 

1.01 

120 

13.0 

6.3 

10.6 

36293 

4101 

397 

1.10*0 

1.40*0 

390.7 

113 

2 

1.36 

1.38 

1.66 

120 

13.0 

6.2 

10.6 

34219 

41834 

277 

1.140*0 

1.00*0 

2603.7 

116 

3 

1.13 

1.60 

1.47 

120 

13.0 

6.2 

10.6 

36229 

41826 

38 

1.170*0 

1.210*0 

2933.9 

115<6> 

6 

1.08 

1.33 

1.62 

120 

13.0 

6.1 

10.6 

34133 

41720 

327 

1.20*0 

1.20*0 

290.8 

120 

1 

1.07 

1.38 

1.43 

100 

13.0 

8.2 

9.9 

30274 

35874 

239 

9.50*0 

9.10*03 

3333.2 

121 

2 

1.99 

2.26 

2.30 

100 

13.0 

8.1 

9.9 

30836 

33949 

244 

9.90*0 

8.80*03 

2733.8 

122 

3 

1.73 

1.97 

2.03 

100 

13.0 

8.1 

9.9 

30247 

33906 

206 

1.020*0 

7.610*03 

230.6 

123 

6 

1.61 

1.87 

1.96 

100 

13.0 

8.1 

9.9 

30Z3S 

3500 

239 

1.060*0 

9.10*03 

201.9 

126 

1 

1.07 

1.46 

1.33 

100 

13.0 

8.1 

9.9 

30263 

33871 

313 

9.30*0 

1.10*0 

390.7 

123 

2 

1.97 

2.23 

2.30 

100 

13.0 

8.1 

9.9 

30263 

33960 

233 

9.90*0 

9.180*03 

280.1 

126 

3 

1.76 

t.96 

2.02 

too 

13.0 

8.1 

9.9 

30829 

33879 

205 

1.020*0 

7.60*0 

2390.6 

127 

4 

1.62 

1.88 

1.93 

100 

13.0 

8.1 

9.9 

30238 

3380 

239 

1.00*0 

9.10*0 

201.9 

43 


entry  Station 
MO.  MO. 

Taaparatura  ot 

Cl)  l  <2)  (3) 

Cdoinot  C) 

Haw 

tot* 

COM) 

Wat if  tulk 
OapC*  T«*> 

<ca>  CC> 

Pr 

m 

HNrautlc  Ofaat 
I*  ta<4> 

ttf 

MU 

Emranca  langth 
(a  Mu 

Mat 

Flux 

<W/«*2) 

12S 

1 

0.31 

O.M 

0.77 

0 

19.2 

8.8 

9.8 

2401 

2707 

2M 

7.491*03 

9.30*0 

1980.7 

129 

2 

1.33 

1.97 

1.M 

0 

19.2 

8.8 

9.8 

2401 

27318 

10 

7.819*03 

7.160*0 

2624.4 

130 

3 

0.99 

1.19 

1.23 

0 

19.2 

8.8 

9.8 

24373 

27310 

10 

7.999*03 

4.910*0 

2607.9 

131 

4 

0.72 

0.97 

1.0 

0 

13.2 

8.8 

9.8 

24373 

2709 

10 

0.10*09 

7.200*0 

278.0 

132 

1 

0.20 

O.M 

0.70 

0 

19.2 

8.8 

9.8 

24373 

2740 

2H 

7.40*03 

1.00*04 

430.4 

1X3 

2 

1.20 

1.43 

1.49 

0 

19.2 

8.9 

9.8 

24322 

274M 

10 

7.790*05 

4.70*0 

S1S.1 

134 

3 

O.M 

1.14 

1.20 

0 

19.2 

8.8 

9.8 

2401 

27478 

iM 

7.900*03 

4.00*0 

2607.9 

133 

4 

0.47 

0.94 

1.01 

0 

13.2 

0.9 

9.8 

24337 

27441 

209 

0.170*03 

7.00*0 

2940.0 

134 

1 

1.73 

2.08 

2.14 

M 

14.5 

10.8 

9.1 

19944 

0990 

223 

4.20*03 

0.010*0 

3763.3 

137 

2 

2.77 

3.04 

3.11 

M 

14.5 

10.7 

9.1 

1990 

0813 

IM 

8.300*03 

7.20*0 

2952.3 

130 

3 

2.94 

2.77 

2.0 

M 

14.9 

10.9 

9.1 

2013 

21883 

199 

8.870*03 

6.00*0 

2499.2 

139 

4 

2.29 

2.94 

2.0 

M 

14.5 

10.9 

9.0 

200(2 

0703 

10 

8.830*09 

7.030*0 

2940.8 

140 

1 

1.04 

1.49 

1.39 

M 

14.5 

10.9 

9.2 

1980 

0338 

230 

8.t40*O3 

9.000*0 

4411.0 

141 

2 

2.00 

2.24 

2.30 

M 

14.9 

10.8 

9.2 

1900 

040 

IM 

8.400*09 

3.900*0 

2624.4 

142 

3 

1.70 

1.89 

1.94 

M 

14.9 

10.9 

9.2 

1970 

03M 

IS 

8.340*09 

4.60*0 

2064.6 

143 

4 

1.31 

1.94 

1.M 

M 

14.5 

10.5 

9.2 

1970 

0314 

143 

8.400*09 

5.410*0 

230.1 

144 

1 

0.48 

1.04 

1.13 

M 

14.3 

10.9 

9.2 

19914 

0309 

20 

8.190*09 

7.30*0 

3873.1 

149 

2 

1.99 

2.24 

2.33 

M 

14.3 

10.9 

9.2 

19934 

219M 

182 

6.490*03 

6.70*0 

202.3 

144 

3 

1.71 

1.91 

1.94 

M 

14.3 

10.3 

9.2 

19934 

0SM 

10 

8.830*09 

4.928*0 

2173.2 

147 

4 

1.09 

1.31 

1.38 

M 

14.3 

10.3 

9.2 

19944 

21337 

10 

6.750*09 

4.10*0 

201.9 

148 

1 

0.40 

0.93 

1.01 

M 

14.3 

10.3 

9.2 

19940 

030 

10 

6.100*09 

6.00*0 

3330.3 

149 

2 

1.99 

2.19 

2.29 

M 

14.3 

10.5 

9.2 

19983 

0374 

181 

6.400*03 

4.00*0 

2624.4 

190 

3 

1.99 

1.74 

1.M 

M 

14.3 

10.9 

9.2 

19911 

0327 

10 

6.610*09 

4.10*0 

1047.2 

191 

4 

1.09 

1.30 

1.34 

M 

14.3 

10.3 

9.2 

1990 

0319 

10 

6.750*09 

5.90*0 

2723.0 

192 

1 

0.92 

0.99 

1.0 

M 

14.3 

10.3 

9.2 

19007 

03M 

293 

6.140*09 

9.10*0 

4626.2 

193 

2 

1.97 

1.82 

1.M 

M 

14.3 

10.3 

9.2 

19813 

21414 

144 

6.430*09 

4.18*0 

Z733.8 

194 

3 

1.09 

1.28 

1.33 

M 

14.3 

10.3 

9.2 

19819 

21300 

118 

6.390*03 

4.40*0 

2064.6 

133 

4 

0.40 

0.84 

0.93 

M 

14.3 

10.2 

9.2 

19784 

030 

m 

4.600*09 

6.00*0 

201.9 

194 

1 

0.31 

0.99 

1.0 

M 

14.3 

10.3 

9.2 

19813 

030 

239 

4.160*05 

9.300*0 

4733.0 

137 

2 

1.99 

1.84 

1.90 

M 

14.3 

10.2 

9.3 

19744 

21332 

187 

4.40*09 

6.20*0 

2733.0 

190 

3 

1.11 

1.29 

1.33 

M 

14.3 

10.2 

9.3 

19781 

0311 

111 

4.340*08 

4.270*0 

1933.9 

199 

4 

0.49 

0.94 

1.0 

M 

M.3 

10.3 

9.2 

197M 

030 

IM 

6.890*09 

5.090*0 

278.0 

140 

1 

1.04 

1.38 

1.44 

72 

14.0 

7.9 

10.1 

22131 

24719 

20 

7.20*0 

1.00*04 

3442.7 

141 

2 

1.87 

2.04 

2.0 

72 

14.0 

7.3 

10.1 

22097 

24730 

IS 

7.300*0 

6.50*0 

1039.0 

142 

3 

1.72 

1.84 

1.0 

72 

14.0 

7.3 

10.1 

2017 

24739 

10 

7.730*0 

9.238*0 

1321.3 

143 

4 

1.31 

1.71 

1.78 

72 

14.0 

7.3 

10.1 

22124 

24739 

IN 

7.930*0 

7.30*0 

210.4 

144 

1 

0.94 

1.17 

1.22 

72 

14.0 

7.8 

10.1 

2210 

24700 

197 

7.230*0 

7.190*0 

2474.3 

149 

2 

1.88 

2.09 

2.10 

72 

M.O 

7.9 

10.1 

2031 

24773 

202 

7.370*0 

7.60*0 

2187.0 

144 

3 

1.44 

1.82 

1.8 7 

72 

14.0 

7.3 

10.1 

2004 

070 

IS 

7.790*0 

4.740*0 

1933.9 

147 

4 

1.94 

1.81 

1.87 

72 

14.0 

7.3 

10.1 

0117 

24737 

243 

7.90*0 

9.60*0 

278.0 

140 

1 

0.48 

0.91 

O.M 

72 

14.0 

7.1 

10.3 

0818 

2420 

20 

7.170*0 

7.460*0 

2474.3 

149 

2 

1.49 

1.M 

1.74 

72 

M.O 

7.0 

10.3 

07M 

2430 

20 

7.470*0 

7.810*0 

2187.0 

170 

3 

1.31 

1.43 

1.49 

72 

M.O 

7.0 

10.3 

0783 

24244 

139 

7.848*0 

9.30*0 

1321.3 

171 

4 

1.10 

1.29 

1.34 

72 

M.O 

8.9 

104 

070 

24204 

187 

7.020*0 

7.40*0 

2069.3 

172 

1 

0.44 

0.90 

O.M 

72 

M.O 

7.0 

10.3 

070 

24208 

217 

7.10*0 

7.90*0 

230.1 

in 

2 

1.31 

1.71 

1.78 

72 

14.0 

7.0 

10.3 

0783 

2420 

211 

7.480*0 

7.910*0 

2107.0 

174 

3 

1.28 

1.44 

1.M 

72 

M.O 

8.9 

104 

070 

2408 

181 

7.630*0 

4.210*0 

1730.4 

179 

4 

1.07 

1.28 

1.33 

72 

M.O 

8.9 

10.3 

0718 

24199 

20 

7.80*0 

8.10*0 

2287.3 

174 

1 

0.13 

O.M 

0.33 

37 

14.0 

8.0 

10.0 

1770 

1001 

28 

9.720*0 

0.10*0 

3335.2 

177 

2 

1.04 

1.24 

1.31 

97 

14.0 

8.2 

9.9 

17890 

19029 

10 

5.900*0 

3.90*0 

218.0 

178 

3 

0.47 

0.81 

0.0 

97 

14.0 

8.2 

9.9 

1780 

18M1 

10 

6.100*0 

4.00*0 

1321.3 

179 

4 

0.29 

0.31 

0.37 

97 

M.O 

8.1 

9.9 

17884 

10073 

140 

6.80*0 

6.30*0 

2396.2 

180 

■0.02 

0.27 

0.34 

37 

14.0 

0.3 

9.9 

17917 

18991 

197 

9.738*0 

7.20*0 

310.0 

181 

2 

0.90 

1.11 

1.14 

97 

M.O 

M 

9.9 

17108 

19018 

IM 

5.970*0 

6.10*0 

22M.4 

182 

3 

0.92 

O.M 

O.M 

37 

M.O 

8.2 

9.9 

178M 

18997 

99 

4.090*0 

3.60*0 

1412.4 

183 

4 

O.M 

0.32 

0.37 

57 

M.O 

8.2 

9.9 

17899 

199M 

18 

4.220*0 

4.90*0 

1960.3 

Him 

bad  t 

■Baraturt  balm  4  d 

•oraaa  C. 

taw  flaw  rata* 

184 

1 

1.44 

1.73 

1.79 

49 

0.9 

0.1 

8.7 

M3M 

(9) 

8 

3.080*0 

9.020*0 

9M.3 

183 

2 

2.12 

2.19 

2.21 

49 

0.9 

0.1 

8.7 

14317 

0 

3.190*0 

7.318*0 

705.3 

184 

3 

2.03 

2.04 

2.07 

49 

0.9 

0.1 

8.7 

1430 

11 

3.278*0 

3.20*0 

30.0 

187 

4 

1.81 

1.M 

1.87 

43 

0.3 

0.1 

8.7 

1430 

IS 

3.350*0 

5.30*0 

344.4 

44 


entry  Station  literature  at 

No.  No.  (O  (2)  <3} 

(dagraea  C) 

Mom 

lata 

<om> 

Uatar  Sulk 
800  T«0 
(cm}  (C) 

fr 

no 

Hydraulic  Olaaatar 
la  te(9)  Hu 

intranet  length 

•a  Nu 

Mat 

Flux 

<W/a*2) 

18t 

1 

0.09 

0.11 

0.12 

0 

27.9 

0.0 

.0 

17999 

12 

3.291*09 

3.90*02 

895.3 

189 

2 

1.11 

1.20 

1.22 

0 

27.9 

33.0 

.0 

17997 

0 

3.381*09 

3.60*0 

909.2 

190 

3 

1.22 

1.29 

1.S 

0 

27.9 

31.0 

.0 

17999 

9 

3.90*09 

1.20*0 

217.3 

191 

9 

0.88 

0.99 

0.97 

0 

27.9 

0.0 

.0 

17939 

19 

3.30*09 

9.371*0 

70.9 

192 

1 

0.19 

0.22 

0.29 

0 

27.9 

0.0 

.0 

17998 

19 

3.20*09 

9.072*0 

03.1 

193 

2 

1.23 

1.33 

1.0 

0 

27.9 

0.0 

.0 

1790 

a 

3.39C*03 

8.20*0 

1093.3 

199 

3 

1.39 

1.38 

1.37 

0 

27.9 

0.0 

.0 

17981 

9 

3.90*05 

1.20*0 

217.3 

199 

9 

0.97 

1.09 

1.0 

0 

27.9 

0.0 

.0 

17981 

19 

3.30*03 

9.30*0 

70.9 

198 

1 

1.98 

1.82 

1.0 

39 

0.3 

0.2 

.3 

12198 

19 

2.20*09 

9.011*0 

893.5 

177 

2 

2.98 

2.87 

2.69 

39 

0.3 

0.2 

.3 

12198 

23 

2.30*09 

8.30*0 

909.2 

190 

3 

2.7B 

2.81 

2.82 

0 

0.3 

0.2 

.3 

12198 

8 

2.40*09 

2.20*0 

328.0 

199 

9 

2.79 

2.0 

2.87 

a 

0.3 

0.2 

.3 

12198 

18 

2.30*05 

9.80*0 

653.5 

200 

1 

0.88 

0.92 

0.93 

10 

29.0 

32.0 

.1 

3708 

9 

8.90*09 

2.991*0 

930.3 

201 

2 

1.88 

1.97 

1.0 

10 

0.0 

32.0 

.1 

3708 

21 

7.20*09 

5.972*0 

909.2 

202 

3 

2.08 

2.10 

2.11 

10 

0.0 

32.0 

.1 

3708 

9 

7.90*09 

1.388*0 

217.3 

203 

9 

2.29 

2.30 

2.32 

10 

0.0 

32.0 

.1 

3708 

19 

7.80*09 

9.20*0 

633.3 

209 

1 

0.83 

0.0 

0.0 

10 

29.9 

32.1 

.1 

3712 

11 

8.911*09 

2.90*0 

516.9 

209 

2 

1.88 

1.97 

1.0 

10 

0.0 

32.1 

.1 

3712 

a 

7.211*09 

3.828*0 

929.5 

206 

3 

2.08 

2.09 

2.10 

10 

0.0 

32.1 

.1 

3712 

9 

7.90*09 

1.221*0 

195.6 

207 

9 

2.22 

2.0 

2.0 

to 

0.0 

32.1 

.1 

3712 

18 

7.80*09 

9.898*0 

718.9 

208 

1 

0.89 

0.0 

0.0 

10 

0.0 

32.2 

.1 

3718 

10 

8.911*09 

2.80*0 

90.9 

209 

2 

1.89 

1.0 

2.0 

10 

0.0 

0.2 

.1 

3718 

21 

7.221*09 

3.998*0 

909.2 

210 

3 

2.09 

2.11 

2.11 

10 

0.0 

32.2 

.1 

3718 

9 

7.938*09 

1.20*0 

206.5 

211 

9 

2.29 

2.31 

2.33 

10 

0.0 

32.2 

.1 

3718 

18 

7.80*09 

9.70*0 

729.8 

212 

1 

1.08 

1.11 

1.12 

a 

0.9 

B.1 

.0 

910 

13 

1.078*08 

3.30*0 

693.3 

213 

2 

2.19 

2.29 

2.0 

a 

0.9 

0.1 

.0 

910 

21 

1.798*0 

3.70*0 

989.2 

219 

3 

2.37 

2.39 

2.90 

a 

0.9 

0.1 

.0 

910 

9 

1.798*0 

1.30*0 

217.3 

219 

9 

2.91 

2.0 

2.0 

a 

0.9 

0.1 

.0 

910 

18 

1.898*0 

9.70*0 

70.9 

APPENDIX  D:  ERROR  ANALYSIS 


NUSSELT  NUMBER 


In  order  to  conduct  an  error  analysis  of  the  Nusselt  number,  an  accuracy  value  for  each  measured 
or  calculated  value  had  to  be  determined.  The  Nusselt  number  was  calculated  from 


Nuh  -  h  Dh 

K 

(Dl) 

where 

Ij  _  ^al(ATalJ 

L  (ATW) 

(D2) 

4  d  (0.3048) 

(2d  +  0.3048  ) 

(D3) 

*AL 

=  conductivity  of  aluminum 

±  13.8  W/m°K 

^AL 

=  temperature  difference  between 
thermistors  in  the  plate 

±0.028°C 

L 

=  distance  between  thermistors 
in  the  plate 

±0.001  m 

ATW 

=  bed  surface  temperature — bulk  water 

temperature 

±0.045°C 

d 

=  water  depth 

±0.001  m 

K 

=  conductivity  of  the  water  at 
the  bulk  temperature 

nil 

REYNOLDS  NUMBER 


Describing  the  Reynolds  number  calculations  similarly, 


ReH  = 


VD  H 
v 


where  V  =  QA 

Q  =  volumetric  flow  rate 
A  =  cross  section  area 
v  =  dynamic  viscosity  at  bulk  temperature 
Dh  =  eq  D3. 


±3.14x  KTWs 

nil 

nil 


Using  the  equation  for  uncertainty  (Holman  1971) 


W  r  = 


1/2 


(D4) 


(D5) 


where  R  is  the  function,  a:, , xr  a'3  ....  xn are  the  independent  variables,  and  W(,  W,,  W3 ....  VV.are  the 
uncertainties  of  the  independent  variables.  The  following  were  calculated  for  two  typical  entries  206 
(low  flow  rate)  and  248  (high  flow  rate). 
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Uncertainty 


Entry  206  Entiy  248 


±10.2 

±157% 

±353 

±15% 

±6.5 

±160% 

±  186 

±15% 

±0.0035 

±50% 

±0.0075 

±  1% 

±1847 

±50% 

±1513 

±  1% 
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APPENDIX  E:  DATA  ANALYSIS  PROGRAM 


'ANALYSIS  OF  FLUME  DATA 
DEFDBL  B,D 

DIM  S(300),T  1  (300),T2(300),F(300),D(300),B(300) 

READ  DATA  FILE 

OPEN  "FDATA"  FOR  INPUT  AS  #1 
INPUT  #1  ,  N 
FOR  1=1  TO  N 

INPUT  #  1 ,  S(I),T1  (I),T2(I),F(I),D(I),B(I) 

NEXT  I 
CLOSE 

500  PRINT  "PRINT  DATA=1" 

PRINT  "ANALYSIS  OF  ALL  DATA  =2" 

PRINT  "RECALCULATE  INDIVIDUAL  DATA  SET  =3" 

PRINT  "EXIT=4" 

INPUT  "SELECT  OPTION";Sl 

IF  SI  =4  THEN  STOP 

IF  SI  =3  THEN  GOTO  660 

IF  SI  =2  THEN  J=1:J1=N:  GOSUB  710 

IFSlol  THEN  GOTO  500 

'*******************************,1,***************************,,, 

590  LPRINT  SPC(  1 3)  "ENTRY  STATION  TEMPERATURE  FLOW  DEPTH  BULK” 
LPRINT  SPC(15)  "NO.  NO.  @.325  @1.59  GPM  CM  TEMP" 

LPRINT 

FOR  1=1  TO  N 

LPRINT  USING  "  ###";I; 

LPRINT  USING  "  #";S(I);, 

LPRINT  USING  "  #.##";Tl(I); 

LPRINT  USING  ”  #.##";T2(I); 

LPRINT  USING  "  ###";F(I); 

LPRINT  USING  "  ##.##";D(I); 

LPRINT  USING  "  ##.##";B(I) 

NEXT  I 

650  GOTO  500 

••m****^*************^**********^*********************************** 

660  J3=0  'RECALCULATES  INDIVIDUAL  OR  GROUPS  OF  DATA 

670  INPUT  "ENTER  FIRST  DATA  SET  NUMBER";  J 
INPUT  "ENTER  LAST  DATA  SET  NUMBER";  J1 
GOSUB  7 10 
GOTO  500 

'*****t*t**************************^********** *********** ***** ******** 

710  'SUBROUTINE  TO  CALCULATE  EXPERIMENTAL  PARAMETERS 
VISC  OF  WATER 

DEF  FNV(T)=H4.375E-06*TA3)+.00 1 3*TA2-.  1 488*T+7.9)*4. 1 339E-04 
COND.  OF  WATER 
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DEFFNK(T)=H  l.ll  E-06*TA2)+.0006889*T+.306)*  1 .7303 

'PR.  NO. 

DEF  FNP(T)=-<  1 .563E-05*TA3)+.004625*TA2-.5 1 78*T+25.6 
'DENSITY  OF  WATER 

DEF  FND2(T)=(999.83952#+ 1 6.945 1 76#*T-.007987040 1  #*TA2-.000046 1 7046 1  #*TA3+ 
.00000010556302 *TA4-2.8054253D-10*TA5)/(  1  +.0 1 687985#*T) 

•ENTRANCE  LENGTH  TO  EACH  STATION 
LI  ( 1  )=  10.55 
L  1(2)=  10.86 
Ll(3)=l  1.164 
L  1(4)=  11.465 

'DISTANCE  BETWEEN  THERMISTORS  AT  EACH  STATION 
L2(l)=.0 12827 
L2(2)=.01262 
L2(3)=.0127 
L2(4)=.01267 

Kl=138  'CONDUCTIVITY  OF  FLUME  BED  W/M  DEG  C 

PRINT  DATA  TO  A  TEXT  FILE 

INPUT  "WRITE  FLOW  CHARACTERISTICS  TO  A  FILE  ";FI$ 

IF  FI$="Y"  OR  FI$=  "Y"  THEN 
CHANGE  THE  NAME  OF  THE  OUTPUT  FILE  HERE 
OPEN  "ADATA.DAT”  FOR  OUTPUT  AS  #2 
END  IF 


INPUT  "REGRESSION  ANALYSIS  (Y/N)";CAL$ 

1220  FOR  I=J  TO  J 1 

IF  I>=1 16  AND  I<=1 19  THEN  GOTO  1300SKIPS  UNSTEADY  DATA 
J3=J3+1 

HF=K1*(T1(I)-T2(I))/L2(S(I))  'HEAT  FLUX 

B 1  =B(I)*  1 .8+32  BULK  TEMP  IN  DEG  F 

SURFACE  TEMPERATURE 

T3=(T2(I)-T1(I))*(0.01905-L2(S(I))-0.003175)/L2(S(I)>+T2(I) 

B2=T3*  1 .8+32  'SURFACE  TEMP  IN  DEG  F 

B3=(Bl+B2)/2  FILM  TEMP  IN  DEG  F 

B9=(B3-32)/l  .8  'FILM  TEMP  IN  DEG  C 


PROPERTIES  BASED  ON  FILM  TEMPERATURE 

B4=FND2(B9)  ’DENSITY  OF  WATER  KG/MA3 

B8=FNV(B3)/B4  'DYNAMIC  VISCOSITY  MA2/S  (FILM  TEMP) 

B5=FNK(B3)  ’CONDUCTIVITY  OF  WATER  W/MA2 

B7=FNP(B3)  PR  NO. 


PROPERTIES  BASED  ON  BULK  TEMPERATURE 

B4B=FND2(B(I))  DENSITY  OF  WATER 

B8B=FNV  (B 1  )/B4B  ’DYNAMIC  VISCOSITY 

B5B=FNK(B  1 )  CONDUCTIVITY  OF  WATER 

B7B=FNP(B  1 )  'PR.  NO. 
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H=HF/(T3-B(I)) 

FM=F(I)*6.30902E-05 

F=FM/{30.4*  D(  1)/ 1 0000) 

H9=4*30.48*D(I)/(2*D(I)+30.48)/10O 

R5=F*H9/B8B 

R6=F*L1(S(I))/B8 

N1=H*H9/B5B 

N2=H*L  1  (S(I))/B5 


HEAT  TRANSFER  COEFF 
'MASS  FLOW  IN  MA3/S 
'FLOW  VELOCITY  IN  M/S 
HYDRAULIC  DIAM  (M) 
RENO.  HYD  DIAM 
'RENO.  ENTRANCE LENGTH 
NU.NO.  HYD  DIAM 
'NUNO.  ENTRANCE  LENGTH 


GOSUB  VELOCITY  'VELOCITY  CORRECTION 

RHV=F*H9/B8B  'RE  .NO.  HYD  DIAM  W/VELOCITY  CORRECTION 


'ADJUST  FOLLOWING  PRINT  STATEMENTS  TO  GET  DESIRED  OUTPUT  FILE 


IF  FI$="Y"  OR  FI$="Y"  THEN 


PRINT  #2,  USING 
PRINT  #2,  USING 
PRINT  #2,  USING 
PRINT  #2,  USING 
PRINT  #2,  USING 
PRINT  #2,  USING 
PRINT  #2.  USING 
PRINT  #2,  USING 
PRINT  #2,  USING 
PRINT  #2,  USING 
PRINT  #2.  USING 
PRINT  #2,  USING 
PRINT  #2,  USING 
PRINT  #2,  USING 
END  IF 


#";S(I); 

#.##";T1(I); 

#.##";T2(I); 

###";F(I); 
##.##";D(I); 
##.##";  B(I); 
###.##";B7B; 

ii  ii ii  it  ii  nit  p  ^ t 

fi  fl  Iftfttrr  y 

######";RHV; 

####";N1; 

###AAAA";R6; 

#.##AAAA";N2; 

TTTrmT .TttT  ,1 II 


###"'!’ 


IF  CAL$="Y"  OR  CAL$="Y"  THEN 
REGRESSION  ANALYSIS 

Y 1  =N  I/(B7BA(  1  /3))  'SETS  Y  OF  REGRESSION  ANALYSIS 

X1=R5  SETS  X  OF  REGRESSION  ANALYSIS 

DM  AT(2, 1  )=LOG(X  1  )+DM  AT(2, 1 ) 

DM  AT(  1 ,2)=DM  AT(2, 1 ) 

DM  AT(2,2)=(LOG(X  1  )A2)+DMAT(2,2) 

RMAT(  1  )=LOG(Y  1  )+RMAT(  I ) 

RMAT(2)=LOG(Y  I  )*LOG(X  1  )+RM  AT(2) 

END  IF 

PRSUM=PRSUM+B7B 


1300  NEXT  I 

DMAT(1,1)=J3 

INPUT  "MORE  DATA  (Y/N)";Y$ 
IF  Y$=”Y"  OR  Y$=  "Y”  THEN 
GOTO  670 

PRINT  “PR  AVG=”;PRSUM/J3 
ELSE 
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IF  CAL$="Y"  OR CAL$=''Y"  THEN  GOSUB  GAUSS 
END  IF 

PRINT  "PR  AVG=";PRSUM/J3  'CALCULATES  AVG  PRANDTL  NO. 

1350  RETURN 

GAUSS:  'SUBROUTINE  GAUSS-SIEDF.L  METHOD  MATRIX  SOLVER 

W=1.8 

XMAT(1)=1  :XMAT(2)=1 

P=1 
SE=0 
1375  M=0 

FOR  1=1  TO  2 
TM=RMAT(I) 

FOR  11=1  TO  2 
IF  1=11  THEN  1400 
TM=TM-XMAT(1 1  )*DM  AT(I,I 1 ) 

1400  NEXT  11 

TM=TM/DMAT(I,I) 

EM=TM-XMAT(I) 

SE=SE+(XMAT(I)/TM)A2 

IF  ABS(EM)<.000001  THEN  1460 

M=M+1 

1460  XMAT(I)=XMAT(I)+W*(TM-XMAT(I)) 

NEXT  I 

IF  P=1  THEN  1500 
CALCULATE  SPECTRAL  RADIUS 
SD=SQR(S1/SE) 

PRINT  "SPECTRAL  RADIUS  ESTIMATE  ";SD 
1500  S1=SE 

IF  M=0  OR  P=2500  THEN  1 580 

SE=0 

P=P+1 

GOTO  1375 

1580  PRINT  EXP(XMAT(1)) 

PRINT  XMAT(2) 

PRINT  "NUMBER  OF  ITERATIONS  REQUIRED  WAS  ";P 
RETURN 


VELOCITY:  'CORRECTION 

NN=.012 
FF=F/.3048 
RT=H9/.3048/4 

FT=8*32. 1 7*(NN/1 .49)^2  *RTA(- 1/3) 

REF=(  F*  H9/B  8  )/FT 

FW=.25264 1 4*REFA-.  1 69 1 2927# 

FB=FT+(2*D(I)/30.48)*(FT-FW) 

UB=SQR(FB*FA2/8) 

REB=(UB*D(I)/100)/B8 

F=UB*(2.5*LOG(3.32*REB)) 

RETURN 

END 


MANNING’S  COEFF. 
'VELOCITY  IN  FT/S 
HYDRAULIC  RADIUS 

'RE(H)/F 
CURVE  FIT 
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